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NOTATION 
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p 
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0 * 

P* 


Axes of the reference triad 

Components of the velocity vector along the 
axes Ox, Oy, Oz 

Absolute temperature 

Potential temperature 

Specific humidity 

Turbulent momentum flux 

Turbulent perceptible enthalpy flux 

Turbulent water vapor flux 

Long wavelength radiation flux 

Latent heat of vaporization of water(2,5 x 10® J,kg“^) 
Karman constant (0.4) 

Acceleration of gravity 
Mass per volume of air 

Mass per volume of atmospheric water vapor 

Specific heat of air ( 

Bowen ratio 

Friction velocity 

Reference temperature 

Specific reference humidity 

Kinematic viscosity coefficient 

Coefficients of molecular dlffuslvlty for 
water vapor and heat 


D and k 



Coefficient of turbulent transfer for momentum, 
heat and water vapor 

Dynamic roughness parameter 

Richardson Number 

Monin Oboukhov length 

Optical thickness 

Modified optical thickness (defined in the 
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THEORETICAL STUDY OP THE INTERACTION BETWEEN 
RADIATIVE TRANSFER AND TURBULENT EXCHANGE IN 
THE LOWER LAYERS OP THE ATMOSPHERE 


Bernard Seguln 


A* 

INTRODUCTION 

Even though It only represents a small part of the earth’s 
atmosphere, the surface layer is of particular Interest. Most 
physical and biological processes which support all animal and 
vegetable life at the surface of the earth in effect depend 
directly on phenomena which characterize the lower layers of the 
atmosphere. Also, disciplines so different such as oceanography, 
agronomy, hydrology, atmospheric pollution studies, etc. are of 
fundamental Interest in this region. In contrast with meteorology, 
which deals with the entire atmosphere, we can speak of micro- 
meteorology for the surface layer. 

This surface layer represents a particular subdivision of the 
planetary boundary layer. As is known, it is the region of the 
atmosphere which is affected by friction caused by the earth's 
surface. It 'is also the region in which the wind direction differs 
from the geostrophlc wind. The surface layer corresponds to the 
first tens of meters of the atmosphere. This region depends 
directly on the adjacent surface below it. The wind direction is 

^Numbers in the margin indicate pagination in foreign text. 
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essentially uniform. It Is useful to consider transforms of 
momentum, heat and mass as being constant along the vertical 
direction. 

It Is sometimes assumed that turbulent phenomena can be con- 
sidered as being statistically stationary for a scale on the 
order of twenty minutes. 

On the micrometeorology scale, the atmospheric motions are 
essentially turbulent and the turbulence has a great Influence on 
all the phenomena. Thus the momentum, heat and mass fluxes, the 
Intensities of which are determined at any time by the energy 
equilibrium established at the surface, are turbulent fluxes. 

The determination of the level of these fluxes is one of the 
main objectives of studies carried on In micrometeorology. 

The direct measurement of these quantities from recording of 
fluctuations of various quantities has the advantage of not 
making use of any hypotheses. Nevertheless, it| encounters experi- 
mental difficulties. Also, usually indirect measurements are 
preferred which use average values of the corresponding quantities. 
Semi-emplrical theories are used and Involve the notion of the 
turbulent transfer coefficient. The accuracy obtained essentially 
depends on the degree to which the various hypotheses used-are. 
valid. These are treated in a special study which make it possible 
to define the experimental conditions which must be maintained In 
order to minimize the errors which are Introduced. 

One of the hypotheses used in general Is that the radiation 
heat flux remains approximately constant through the surface 
layer and consequently does not Interfere with the turbulent heat 
flux which can then be calculated Independently. Even though 
this hypothesis Is explicitly or Implicitly assumed in most work 
on the turbulent structure of the surface layer. It nevertheless 
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has been seriously placed /ln doubt In recent studies. 


The work which we will present here attempts to estimate the 
variations in radiation transfer through the surface layer of the 
atmosphere. We will estimate their effects on the distribution 
of turbulent heat fluxes in this region. 

The present paper has the following plan: 

- After recapitulating various simplifying hypotheses which 
lead to assuming that the vertical turbulent fluxes are 
conserved (Chapter I), the method of determining these fluxe 
by seml-emplrlcal theories is treated (Chapter II). We 
will stress the hypothesis which Involves the radiation 
transfer, andjthe steps to control it are presented: the 

radiation transfer is Initially assumed to be constant and 
then its vertical distribution is calculated from theoreti- 
cal temperature and humidity profiles which result. This 
makes it possible to determine whether the initial hypothe- 
sis is valid or not (Chapter III) . The theoretical 
temperature and humidity profiles are studied first of 
all (Chapter IV). Then the method of determining the 
divergence of the radiation flux from these profiles is 
presented (Chapter V) . The details of the numerical cal- 
culations as well as their results are given (Chapter VI). 

An analytic method makes it possible to rapidly estimate the 
radiation transfer effects (Chapter VII). The general 
conclusions are finally given in (Chapter VIII).* 


*Thls paper was done within the more general framework of 
atmospheric-ocean interaction studies carried out at the I.M.S.T. 
Only the particular case of the ocean was considered in numerical 
applications. The results would be applicable to an arbitrary 
natural surface. 
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I. THE CLASSICAL THEORY AND NOTION OF ” CONSTANT FLUX LAYER” /Jl 

The theory of the turbulent surface layer of the atmosphere 
has been discussed In nuj^erous works of which we would like to 
mention Sutton (1953) j Priestley (1959) j Lumley and Panofsky (1964), 
Monln and Yaglom (1966), Munn (1966), Zllltlnkevltch (1970). 

In the particular case of the ocean-atmosphere Interface, 
the basic equations are presented and discussed In the article of 
Coantlc ( 1969 ) which Inspired the following paragraphs. 

II . Simplified equations applicable to the atmosphere 
surface layer 

The averaged equations applicable to the surface layer are 
established using the general statistical equations of atmosphere 
mechanics. The particular properties of the region under consi- 
deration will be taken Into account: 

- The variations and fluctuations of the specific mass are 
sufficiently small so that they do not Influence either 
the definition of the macroscopic quantities (Havre 1965) 
nor the form of the averaged equations which follow. 

- The effects of Coriolis forces can be assumed to be negli- 
gible. The Interaction phenomena are then two dimensional 
on the average. 

- These phenomena have a boundary layer character: the 

variation of the various variables Is much more rapid along 
the vertical direction than along the direction of motion. 
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The following simplifications result: 


V =. O 

^( 5 ^) » 


( 1 ) 

( 2 ) 

(3) 


The equations which control the turbulent transfers of 
momentum, water vapor and enthalpy in the lower layers of the 
atmosphere are then written as: 


llh'* 

L -^x Si 




dz: 




^ + U i + w ^ 

^z. 


- K' + Rl 


(H) 


(5) 


( 6 ) 


/4 


and the equation of continuity takes on the simple form 


•^u 


c>Z 


(7) 


The simplifications Introduced up to the present are generally 
valid. However, the system consisting of equations (4), (5), (6) 
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and (7) is too complex to be easily solved. It can be simplified 
using additional hypotheses and of course there will be loss of 
generality. 


12. Hypothesis of horlziontal homogeneity 


The processes studied are assumed to occur over horizontal 
ranges which are sufficiently large so that the phenomena can be 
considered to be "completely developed". In other words the 
averaged derivatives along the direction Ox can be completely 
Ignored 


SX 


(x)=o 


( 8 ) 


In the case of the earth's surface, this hypothesis is some- 
times open to question because of the existence of various 
heterogeneous features of various types. Numerous studies, 
theoretical as well as experimental, have been carried out recent- 
ly on the subject. They make it possible to delimit the range of 
validity of the preceding hypothesis or the minimal horizontal 
distance in the wind direction from a heterogeneous region and 
for a given observation height (see for example Elliott 1958, 

Dyer 1963, Rider, Philip and Bradley 1963, Panofsky and Townsend 
1964, Peterson 1969 etc...). In the particular case of the 
ocean, if we exclude coastal regions and contact regions of 
marine currents having different thermal properties, this hypothe- 
sis can be used much more extensively, because of the homogeneous 
nature of the ocean surface. ^ 

The advectlon terms in the first term of the preceding 
equations can then be Ignored and the system is reduced to 
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Tt ) = - ^ {-'"^ * ^] 

13. V steady state hypothesis 

Even though the question Is still being discussed today, it 
is generally assumed that the spectrum of dynamic turbulence in 
the atmospheric surface layer is separated into two distinct parts 
by a spectral gap which extends from a few tens of minutes to 
several hours (see, forj example, Lumlej^and Panofsky 1964 , Panofsky 
1969). These phenomena can consequently be considered as statis- 
tically steady over a time Interval which varies from five to 
thirty minutes, depending on atmospheric conditions. Also, 
except for transition periods (such as during the rising and 
setting of the sun or when clouds pass by) the fundamental calori- 
metric input produced by the sun’s radiation varies little, and 
the energy equilibrium can be considered as being established. 

We can therefore write: 



( 9 ) 


( 10 ) 


( 11 ) 
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and the system of equations becomes 


^ 1- R.’5h+ puV I = "^=0 
■&*. I ^ az. ' J ■&Z- 


( 13 ) 


l-z 


( 14 ) 




|^{s-R] = o 


( 15 ) 


/± 

The total momentum flux t, the water vapor flux J and the 
enthalpy flux S + R , are then constant along the vertical 
direction. 


Over almost all of the surface layer the transfer by molecular 
diffusion can be Ignored with respect to the turbulent transport, 
and we can therefore write 


L = ^ u V = to 
T ^ f — 'Jo 


( 16 ) 

( 17 ) 


S + R 


Sq + Rq 


with 



(18) 
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l4 . Hypothesis regarding the conservation of radiation flux 


The validity of the hypotheses which we have mentioned can 
be checked In each particular case, as we have seen. For example 
it Is sufficient to verify that the selected space and time regions 
have been appropriately selected. 

The system of equations (l6), (17), (l8) Is then strictly 
valid. There are no turbulent quantities Involved except for the 
radiation flux R, which has a different fundamental nature and 
cannot be determined using the same methods. In a general way, 
more or less Implicitly, this radiation flux Is assumed to be 
constant over the thickness of the surface layer. This hypothesis, 
which Is justified because of the small atmospheric thickness 
(a few tens of meters) considered, as well as the high transparency 
of humid air, seems to be corroborated by the agreement between 
the theoretical positions and the experimental results which are 
most often observed. 

It follows that the turbule'nt fluxes of momentum, water vapor, 
and enthalpy are In general considered as constants along the 
vertical direction 


r = P. 

T = 

S s 5 

o 


( 16 ) 


(17) 

(18) 
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This Is why the surface layer Is often called the "constant 
flux layer" . 

II. SEMI-EMPIRICAL RELATIONSHIPS BETWEEN THE PROFILES AND 
THE FLUXES 

21 . Methods of determining the turbulent fluxes 

There are numerous- such methods which use various techniques 
(see Readings, etfal ,|1969) . Only two of them use the properties 
of turbulent flows. We will not mention the others here without 
evaluating their correctness. 

In a general way. If is any one of the properties related 
to mass, the turbulent flux corresponding to It Is of the form 

4[i) . ^ 7i' tw) 


Its value can-' thus be obtained by simultaneously recording 
instantaneous fluctuations of the property pc\ and the vertical 
component w at an arbitrary height z. The problem Is essentially 
an experimental one connected with accuracy, response time and 
space definition of the transducers. The measurement of Reynolds 
stresses by this procedure Is a classical one. The measurement 
of the heat flux is relatively well known (Buslnger et al.jl967). 

On the other hand, the determination of the water vapor flux 
(such as the carbon dioxide flux) is a delicate problem because 
of the difficulty of measuring humidity fluctuations (Dyer and 
Maher 1965, Miyake and Mac Bean 1970, Leducq 1970). 

In spite of progress made in this direction, the required 
instruments are complex and expensive and are relatively delicate,| 
which limits their use at the present time. This explains the 
10 



development of other methods which require the Introduction of 
additional hypotheses but which can be used In a relatively 
simple way. The fluxes are simply related] to the vertical 
profiles of the average values of the corresponding properties 
by means of the turbulent transfer coefficients. These coeffi- 
cients, defined In analogy to the molecule transfer coefficients, 
are determined In part by theoretical considerations and In part 
by empirical extensions which are the results of experimental 
data. This explains the name seml-emplrlcal given to the method 
In question. 

22 . Expression for the turbulent transfer coefficients 

The various coefficients are defined as /8 




T = - p ICe 

' "Sz. 


S =- 




( 20 ) 


( 21 ) 


( 22 ) 


T|he averaged values are obtained over conveniently selected time 
Intervals (see paragraph 13). 

kj^, kjj and k^ are turbulent transfer coefficients for the 
momentum, heat and water vapor, respectively, and we will' briefly 
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recall the presently used expressions. 


In the first approximation, we consider the case of adiabatic 
conditions. The gravity forces, or If we prefer, the effects of 
thermal stratification, are then negligible (in the surface layer, 
the adiabatic conditions are essentially equivalent to the thermal 
neutrality conditions defined by * 0) . 

The situation Is analogous to the situation encountered In the 
laboratory In the internal part of the boundary layer over a 
rough plate. A logarithmic velocity distribution in the vertical 
direction results which Is Identical to the classical "wall law". 

and, by integration (23) 

( 24 ) 


J,21 

u (z) =. Lcci ^ 


is an Integration constant called the roughness parameter 
which translates the overall dynamic effects which are the result 
of the nature, height and spatial distribution of the surface 
uneveness . 


Under adiabatic conditions, the various turbulent transfer 
processes depend only on the dynamic process: the various 

properties are transported In a passive way. In general It Is 
assumed that the transfer coefficients are related by the simple 
relationships 






H 


and 






*Translator ' s note: Figures missing in foreign text. 
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and 


are the Inverses of the "turbulent Prandtl and /9. 


O' 

where 
Schmidt numbers". 

In analogy with the friction velocity Ujj , we can define a 
reference temperature and a reference humidity Pjf, such that 








S 




(25) 


(26) 


The temperature and humidity profiles are then written in the 
form (Zllltinkevltch 1970, Monin 1970) 




6o = -^ La,i-+S©, 

Zo I 




t" o — i— oq -t- ^ "fo 


(27) 


( 28 ) 


the quantities ^eJ, and .s'y translate the deviation which exists 
between the mechanical energy transfer mechanisms, ~6n| the one hand,] 
and those of other turbulent transfers , dh| the other hand7”near| the 
surface when it is very rough. Stated differently, the roughness 
parameter z^, which is defined with a dynamic process, is not] ' 
necessarily identical to those which could be defined for other 
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transfer processes. 


In effect, in the surface layer, the thermal neutrality 
conditions are rare. The temperature gradient Is In general 
over-adlabatlc during the day and there Is Inversion at night. 
Also, In the general case, the gravity force effects affect the 
structure of turbulent flow. The now classical works of Monln 
and Obukhov (again published by Monin and Yaglom 1966) made it 
possible to perform a strict analysis of the surface layer In the 
general case. Without going Into details, we would like to state 
that these authors defined three external parameters In which any 
Internal variable of the surface layer can be determined: these 

are the quantities Uj^, and L. 


L Is the Monln-Obukhov length. The expression for it Is the 
following and Is established by similitude laws 




(29) 


The non-dimensional ratio is a variable which is /20 

characteristic for the state of stability of the atmosphere at 
a given level. It is substituted for the more classical Richardson 
number, Rl, to which it can be related. 


We can write 


l4 



(30) 

(31) 

(32) 



This system Is equivalent to the following system for the 
turbulent transfer coefficients 







_ 2L 

(^) 


( 33 ) 

( 34 ) 

( 35 ) 


In the particular case of adiabatic conditions, which we 
considered previously, and — >0|. 

This makes it possible to determine the following conditions 
for the functions <j> m» 4^H|and <|De' 


(o) = 
(o) = 

4>r I ®) = 


d 


K 


( 36 ) 

( 37 ) 

( 38 ) 


The functions c|> |v|, h| c|5E|could only be determined /II 

theoretically up to the present. 
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Their variations as a function of the parameter^' 
established from experimental results. 

By adopting a seml-emplrlcal relationship for each of the 
functions and It becomes possible to directly 

determine the fluxes from the vertical profiles of the correspond- 
ing averaged quantities. This method Is now used. Except for 
experimental problems. Its accuracy depends essentially on the 
validity of the seml-emplrlcal relationships used. Therefore 
it can be understood why the possible vertical variation of the 
radiation flux Is of Interest, not only during measurement 
taking but also when the mlcrometeorologlcal measurements are 
Interpreted . 

Note : An extension to the definition of the Monin-Obukhov 

length must be made when the evaporation flux Is large. In order 
to take Into account the gravity effects which are the result 
of humidity gradients. The effects of stratification only appear 
as correction terms as we will see In the following. We did not 
take Into account this generalization. /12 

III. PROBLEMS CAUSED BY THE VARIATION OF THE RADIATION FLUX 
AND THE METHOD PROPOSED FOR AVOIDING THEM. 

31 . General considerations on the radiation transfer In the 
atmospheric surface layer 


'Can only b^ 


As we have seen, the hypothesis of conservation of the verti- 
cal radiation flux Is a basic hypothesis of the classical theory 
of turbulent transfers In the vicinity of the surface. The lower 
layers of the atmosphere constitute an absorbing medium for long 
wavelengths, because of the presence of relatively large concen- 
trations of water vapor and carbon dioxide. The presence of 
temperature gradients brings about a certain variation In the long 
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wavelength radiation flux. The determination of such effects 
is done in a classical way in meteorology (atmosphere thickness 
on the order of 500 m) . A sufficient accuracy has been obtained 
by using diagrams (Elsasser, Kew, Yamamoto, etc...). ]^t is mor^ 

delicate in the region of the surface layer because of the small 
thicknesses considered. We were able to collect a certain number 
of papers on this question and we would like to mention the follow- 
ing: Brooks D. L. (1950), Deacon (1950), Robinson (1950), Yama- 

moto and Kondo (1959) j Punk (I960), Punk (196I), Lalkhimann et al 
(1961), Gaevskaya et al (1963), Elliott (1964), Godson (1965), 
Hamilton (1965), Kondratiev^ (1965) , Adunkowskl and Johnson (1965), 
Atwater (1966), Lleske and Stroscheln (1966), Seo et al.J(1968), 
Hlnzpeter and Heinrich (1969), Paraponova (I969). 

Based on these numerous papers, one could think that this 
problem had already been sufficiently studied. In reality, the 
large majority of the studies mentioned approached the question 
from a different angle, because an attempt is made to predict the 
variation of a non-steady thermal regime from certain experimental 
data . 


In effect, within the surface layer assumed to be homogeneous 
in the horizontal direction, the equation which governs the 
turbulent enthalpy transfer can be written as 


where 


S • -VS* tt 1 “ 


(39) 


In the general case of non-steady conditions, the heating 
and cooling rate of the lower layers results in the combination 
of the divergences of the radiation and turbulent fluxes. 
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As far as the nocturnal cooling is concerned, the prediction /13_ 
of which is extremely Important for economic reasons (prevention 
of ice, fog, etc.) and which is in general associated with small 
turbulent convections, can be considered negligible with 


respect to 


■^21 


We then have 




(40) 


The nocturnal cooling rate is therefore directly proportional 
to the radiation flux divergence and this is one of the reasons 
for the relatively large number of papers devoted to this question. 


For the non-steady regime, the vertical variation of the 

turbulent flux can be derived experimentally from the equation 

(39), if M] Is determined and is measured. Robinson (1950), 

Funk (i 960 ), Elliott (1964) also compared the true cooling 

yl ^ 


to the radiation cooling 
to the term -^^1 




The difference corresponds 


We can therefore see that there are various ways of conslder- 
jng the problem. For our part, we proceed as follows: our goal 

as to verify the validity of the hypothesis made regarding the 
radiation transfer within the framework of the classical method 
of determining the turbulent fluxes. We are| therefore led to 


assume that a steady regime exists. 


property that 


as 

>z. 


~2>z; 


will then have the 


(41) 


It is then clear that the method of determining the divergence 
of the radiation flux can be used in the non-steady regime. 
Consequently the numerical values obtained can be compared with 
18 



those obtained by various authors which we mentioned. 

32 . Influence of the variation of the long wavelength 
radiation flux on the quantities defined In the 
constant flux surface layer 


For a long time. In the study of turbulent transfer near 
the surface. It has been assumed that the divergence of the 
radiation flux could be considered negligible because of the 
small atmospheric thickness under consideration. This opinion 
was confirmed by the conclusions of Robinson (1950) who, according 
to experimental results, was led to believe that was negligible 

except possibly below 1 m. Nevertheless, wg' should note that the 
quantity considered In general was the relative variation of the 
radiation flux proper. This quantity Is effectively very small. 
However, for the problem with which we are dealing, this Is not 
the quantity we are concerned about. Instead we are concerned 
with the divergence of the radiation flux relative to the turbu- 
lent heat flux. /l4 


In effect equation (17) can be written as 


where S (z') = So with R (z) - R© A 


(^ 2 ) 


It Is therefore clear that the turbulent flux 5(2)! cannot 
be considered as a constant and equal to except to the extent 
that Is negligible with respect to S^. A small rela- 
tive variation In the radiation flux can lead to a 

value of which Is not negligible with respect to 1 If 

Is small with respect to . It Is therefore necessary to simul- 
taneously know the quantities which characterize the radiation 
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transfer on the one hand^ and) the turbulent transfer on the other 
hand. This justifies the Interaction term which we have used in 
this section. 


Under the conditions under discussion, the variation in the 

radiation flux R brings about a variation in the turbulent flux 

'S, which then is assumed to be constant and equal to S.- . Thus 

o 

an erroneous estimation of the flux P results when it is determined,] 

if the temperature gradient and the coefficient k„ are known. 

n 

Inversely, an erroneous estimation of kj^ results when it is derived 
from measurements of the flux and of the gradient. In effect, 
strictly speaking kj^ is defined 


S{z) 




k I 


H 


■2>Z. 


( 43 ) 


If S is assumed to be constant, in reality we determine a 
global coefficient k„ such that 








( 44 ) 


We can see that the relative error of k„ is equal to that 

n 

of S. 



The ratio is used j'ointly for determining the relation- 
ship > which Introduces an error due to the influence of 

radiation transfer. 
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A certain number of authors have for a long time Indicated /I 5 
the possibility of radiation transfer (Priestley 1959, Mac Vehll 
1964, Lumley and Panofsky 1964, Godson 1965, Webb 1970, Oke 1970). 

Certain authors based their calculations on experimentally 
determined profiles and detected this effect In certain particular 
mlcrometeorologlcal situations. 

For the case which he studied, Elliott (1964) calculated a 
variation of l4^ In the turbulent flux, even though the wind 
velocity was quite high (7 to 8 m/s). He also adds: 

"However, much stronger Inversions than the ones on which these 
computations are based are frequently found. These stronger 
Inversions would lead to larger radiative cooling rates and 
quite likely to eddy flux divergences which could not be Ig- 
nored In considering such low-level phenomena as fog-forma- 
tion" . 

Also, Kraus, referenced by Munn (1966), studied three parti- 
cular cases for which he found a variation of 2 %, 7 % and 50^, 
respectively . 

All these observations led Munn (1966) to write: 

"The observed temperature structure In the boundary layer Is 
the Integrated result of turbulent mixing and radiative trans- 
fer'.^ One of the unsolved problems of micrometeorology Is the 
determination of the relative Importance of the two components". 
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33. Description of the method utilized for quantitatively /l^ 
estimating the Influence of radiation transfer 


The Influence of radiation transfer seems to not have yet 
been studied In the most general case, where the various parameters 
are not determined numerically In advance and can take on arbitrary 
values corresponding to any mlcrometeorologlcal situation. This 
global study has two advantages: 

- on the one hand, the Influence of the various parameters 
can be shown analytically, which makes It possible to 
better evaluate the relative Importance of the various 
processes . 

- on the other hand, generally valid equations could be then 
applied to each particular case. For each determination 
of the flux or of the turbulent transfer coefficient, 

an approximate evaluation of the error Introduced by 
Ignoring the radiation transfer can be obtained relatively 
easily. This evaluation seems to be necessary at the 
present time, as the studies on turbulent transfer 
mode recently by Webb (1970) and 6ke (1970) have shown.] 

These two authors deplore the absence or Inaccuracy of 
methods which would make It possible to evaluate the Inf lu^nc^e] 
of the radiation transfer during their experiments. 

We are therefore forced to carry out a general study. It Is 
based on the following scheme which was conceived by M. Coantlc 
and which we personally developed while collaborating with him. 

As generally assumed (and this Is the point which we want to 
verify), the surface layer Is first assumed to be a "constant 
flux layer" and the radiation flux variation Is assumed to be 

negligible. According to hypothesis, the turbulent flux S Is then 
constant . 
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It Is then possible to establish theoretical temperature and 
humidity profiles. If the fundamental parameters are given from 
the beginning and If the usual seml-emplrlcal relationships are 
adopted. This method of calculation, which Is a completely 
classical one. Is only valid If S Is effectively constant. Prom 
the temperature and humidity profiles It Is possible to calculate 
the divergence of the radiation flux. By Integration Its variation 
over the surface can be calculated. It Is assumed that steady 
conditions prevail. The total enthalpy flux must then remain 
effectively constant and the absolute value of the variation 
obtained Is equal to the variation of the turbulent flux S. Thus 
S was assumed to be constant. If the variation given by the 
calculation Is small, l.e. on the order of magnitude of the errors 
usually made In micrometeorology (on the order of 10^), the Initial 
hypothesis will remain approximately valid. On the other hand. 

If the variation which results Is higher, the hypothesis Is In- 
validated by the conclusions ;| the error Introduced by removing 
the Influence of the radiation transfer cannot be considered as 
negligible. The temperature and humidity profiles are different 
from the ones which we used. The general solution of the problem /17 
then requires recourse to a method of succeslve approximations. 

The study then consists of two main parts: 

1. Establishment of theoretical temperature and humidity 
profiles In a constant flux layer 

2. Calculation of the divergence and variation of the long 
wavelength radiation flux. 
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IV. DETERMINATION OP THE VERTICAL DISTRIBUTION OF TEMPERATURE /l8 
AND HUMIDITY FOR THE HYPOTHESIS OF A CONSTANT FLUX 
LAYER 

4l . Parameters and relationships required for defining the 
profiles 

The various transfer modes all depend directly on the dynamic 
process. The determination of the velocity profiles is Indlspen- 
slble for defining the vertical temperature and humidity profiles. 

In the constant flux layer It Is known that the dynamic 
process is defined by specifying three fundamental parameters 
Ujf, and L (Lumley and Panofsky 1964). 

In the particular case of the ocean, and for the problem 
which we are treating here, It is generally assumed that the 
roughness of the Interface Is determined by the value of the 
friction velocity u^ (the statistical characteristics of the 
waves are assumed to be determined by the local wind velocity 
as the first approximation) . 


z Is then a function of u«' which we will now examine. The 
o * 

number of fundamental parameters Is reduced to two: u* and L. 


4ll . The method of calculating the velocity profile Is then 
as follows :| ' ■ ■ 


plrom ( 3 ) we have 



(45) 


by Integration, 







(46) 
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At the Interface, by definition of Zg 

2 > 2 J . 


u(z) — >0 when 


Thus , 








(47) 


The velocity profile is thus completely determined by 
specifying the pair u^, L if the functions (u*) and 
are known. 

4l2 . As far as the thermal process is concerned, equation /19 
( 31 ) makes it possible to wrltej 

: e, 

“^>2: Z 


+" 


(48) 


or by integration 


e(z)= e* Jiil%)dz.+ c 




The Interface equation can be written as (see paragraph 22) 
© when Z~ 


which leads to the expression for the temperature profile 


6(zl) - 0o - r -^clz 


(49) 


It is therefore determined from the velocity profile by 
specifying _(^_§0o and 0* . In reality, the parameter 

is already known. In effect, is defined by the relation- 
ship (26) According to the definition of the 

Monln-Obukhov length • Except for very pronounced 
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advection conditions, which are then eliminated according to 
hypothesis in our study, the natural exchange mechanism only 
results in moderate temperature gradients in the surface layer. 
0 can be replaced by 0^ as an approximation to within \% . We 
then have 

s _ _ u^l^ Go 


By combining (26) and (48), we can then obtain 0^: 
0^ is then determined if u^, L and 0^ are known. 


(51) 


Thus, for the thermal process, the velocity profile is 
assumed to be known and the determination of the temperature 


profile requires that the quantities 
Ht)\ be specified. 




and the function 


413 • Finally, the evaportionl-process is defined in a 
analogous way. According to (32) in effect we have 


/20 


or 



(52) 

(53) 


with the condition at the Interface t 




Srlwhen 
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But F* Is simply related to by the Bowen ratio 


3 = -§- 


In effect, by definition 


or 





T 

fhnyr 

S 

Z_ 

5 ' 


which results In 


u^Bo 

J^3 


( 54 ) 


(25) 

( 26 ) 


(55) 


The humidity profile Is therefore determined If the velocity 
and temperature profiles are known by specifying 

and the Bowen ratio. 

In conclusion, the various profiles can be completely 
calculated In any case If the quantities L, 0 ^f’ol 

and B are specified and If the seml-emplrlcal relationships 

and are known. We will now examine 

these various details. 

42. The roughness parameters on the surface of the ocean /21 

Because the dynamic process Is predominant. It Is necessary 
first of all to define the dynamic roughness parameter . This 
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parameter, which Is an Integration constant which results from 
extrapolating the velocity profile. Is a function of shape, height 
and spacing of ob stacles wh ich make up the roughness In the ca se| 
where the surface Is fixed. In the case of the ocean, the surface 
Is deformable under the action of the wind, which leads to the 
development of waves. The waves have variable shapes and speeds 
In time and space. This brings about extremely complex dynamic 
effects which are not well known. In order to apply semi-emplrl- 
cal theories of atmospheric turbulence. It Is generally assumed 
that these effects can be translated by means of the parameter 

which only depends on the local dynamic interaction, represented 
by the friction velocity u^^. All the information for this rough- 
ness parameter ls{ of an experimental nature. The first experi- 
ments were carried out In the ocean, but It seems very difficult 
to Interpret the results, as Figure 1 shows, taken from the 
book of Roll ( 1965 ) . 

The research carried out In the wind tunnel (Hldy and Plate 
1967 j Shemdln I 967 , Karakl and Hsu 1968, Wu 1968 and 1969 j etc...)] 
made It possible to better understand this phenomenon. In the 
development of the various waves which Influence the dynamic 
roughness of the ocean surface, there] is not only the local tur- 
bulent stress caused by the wind and translated by the quantity 
Ujf, but there Is also the "fetch" X ("fetch" Is a nautical term 
which designates the surface traversed by the wind'. . Thus In the 
case of the ocean, the distance to the shore or the origin of 
the perturbation). For the same wind velocity, the wind dimentlon 
Increases with X (Figure 2). When X Is small, u^^ also varies as 
a function of Xi”be~causejof the proximity of the "attack band" 
there is a complex variation of z^ as a function of X (Figure 3). 

In the present study, the flow Is assumed to be completely 
developed and the role of "fetch" can be completely eliminated 
according to hypothesis. It will therefore be necessary. In order 
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to extend the results obtained in the wind tunnel to the problem 
which we are considering, to isolate the only variation of as 
a function of u^f for a fetch assumed to be Infinite (in fact, a 
fetch on the order of 10 km seems to be sufficient) . 


For X given^ two different flow regimes are observed 
corresponding to the wave dimensions which make up the roughness. 
For small u*, the average heights of roughness are small and the 
flow regime corresponds to the "hydraulically smooth" regime, 

Is then connected with the thickness of the viscous sublayer and 


has the form 


Zo- Q 


(Monin 1970). For higher u« values, the 

“ 1r 

regime tends to become "hydraulically rough" after a transition 
region. Charnock established a theoretical relationship of the 


form . b depends on the value of the fetch according 

to the previous developments. 


Overall this scheme has been confirmed by experimental /22 

results Hidy and Plate and Shemdln (Figures 4 and 5) . 


We should like to state that these authors obtained a limited 
"fetch" on the order of magnitude of the wind tunnel length. 

The extension to a practically Infinite "fetch" corresponding 
to the case of completely developed flow over the surface of the 
ocean was done by Wu (1969) (Figures 6 and 7). The results ob- 
tained are In agreement with the qualitative model which resulted 
from the wind tunnel studies and were the results of various 
observations.* 


*Thls agreement is different from the dispersion found by 
Roll. First of all the results collected by Roll are much older 
than those considered by Wu. Also the measurement accuracy has 
been considerably increased. The theories for the velocity pro- 
file shapes has advanced. In particular, the error associated with 
thermal stratification has been eliminated. (Continued on next page) 
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We therefore utilized the numerical values of from 
Figure 7. 

r~ \ 

For thermal and mass exchange processes the terms and 

translate the deviation which occurs close to the rough 
surface between the mechanical energy and the various other trans- 
fers. Various recent studies (Barry 1965, Chamberlain 1968, Cowan 
1968, Thom 1968, Zllltlnkevltch 1970) have shown that the terms 
could be correlated with the ’’roughness Reynolds number” *1 

and §foj are Increasing functions of the ’’macrovlscoslty ” 
u^ z^ defined by Sutton (1953)* Consequently they must not be 
considered except for very rough surfaces, as Chamberlain (1968) 
showed In experiments. 

The small 'va^luesj of the roughness parameter z^ for the ocean 
surface then leads to negligible values of and S'fJ and we 

therefore write 


= o 


/23 

43 . Determination of the surface temperature and humidity 
and the Bowen ratio. 


It would first seem Illogical to regroup the surface humidity 
and the Bowen ratio. These quantities have a common feature 
which Is Important for a study: strictly speaking, they cannot 

be considered as given quantities of the problem but Instead 


*( Continued from preceding page) Finally, the role of the "fetch”, 
which could be great for measurements near the coasts, has been 
considered as we showed. These various reasons .lead us to believe 
that the measurement dispersions observed by Roll are due not to 
the phenomenon Itself but to Inaccuracies In the experimental 
determinations of z and to the very different fetch effects. 
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must be considered as results. 

As we stressed In the Introduction, the distribution of 
energy between the various fluxes results In an equilibrium which 
Is established at any time at the level of the surface because 
of the action of the heat Input from the solar radiation. 

Depending on the boundary conditions, in the air and also 
in the water, (wind velocity, air, temperature and humidity above 
the Internal boundary layer, water temperature with depth) and 
the solar radiation absorption mechanism In the water, the system 
leads to a unique solution and the quantities u*, L, ©tjoi fo] and 
B as well as S and J are the results. The rigorous solution of 
the problem requires the complete solution of the system. We will 
consider the first part as being solved and will adopt plausible 
values of these parameters. 

Among these, three are independently specified: Uj^, L and] 

®o- 

We will assume that F^, the surface humidity. Is the saturating 
humidity of sweet water at temperature 0^ and Is therefore com- 
pletely determined by specifying this parameter. F^ also depends 
on the radius of curvature of the surface and the salinity of the 
water; also] a surface film and molecular scale effect can have an 
Influence. However, In the first approximation, the Influence 
of these factors can be considered to be negligible. 

The Bowen ratio B Is more difficult to specify. It Is not 
completely Independent of the other parameters. Below the ocean 
Its sign Is generally determined by the sign of the heat flux S. 

The flux J Is almost always directed upward (evaporation takes 
place). J Is therefore positive except In an Inversion situation. 

S Is '<1 o] In an inversion situation. S is > oj under 
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over-adlabatlc conditions. Thus B = 

•b 1 

the sign opposite to that of L. 


has the sign of S or 


We could then consider specifying various arbitrary values 
of B, which would be positive for L < Dj and negative for b > o| 

B would then be an additional variable, comparable with u^^, L and 
0Q. But the number of possible combinations between these various 
variables then becomes very large. We preferred utilizing the 
values obtained by Gordon and reported by Roll (I 965 ) (page 25^) 
where B was experimentally correlated with the temperature 
difference between the air and the surface. Even though there Is 
no experimental justification for this correlation, we made use of 
It In this study In order to simplify our analysis. /24 


44. Discussion of th e seml-emplrlcal relationships 
and 


This section alone could_be the_topic of a paper which]' 
would be as Important as this entire paper. But this Is not our 
purpose and we will give an outline of the problem and specify 
the expressions which we use. 


First of all It should be stated that, by] adapting a seml- 
emplrlcal relationship of any kind, this does not fundamentally 
modify the significance of our results. As we will see In the 
following, to first order, the stratification effect Is negligible 
and the adiabatic approximation H ^ | is sufficient.] 

The form of the seml-emplrlcal relationships does not enter except 
as correction factors. 


The seml-emplrlcal relationships ^ h| and 4 ^e: 1 translate 

the Influence of the thermal stratification to mechanical energy 
transfers, heat transfer, and mass transfer, respectively. The 
use of similitude relationships inalces| use of the Monln-Obukhov 
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length L. It became possible to replace the empirical relation- 
ships formally based on the Richardson number R1 [Jsee for example 
Deacon (19^8)3 by more satisfactory theoretical expressions. 

Since the study of Monln and Obukhov date from 195^ s these are 
rather recent papers. It seems that a coherent theory Is develop- 
ing and we will now outline it. 


44l . We will start by studying the relationship 
fined by the equation (30) 


- 4. z. 




de- 


Por adiabatic conditions (see paragraph 22), l.e. 

~ — >0 — >l| . Monln and Obukhov considered the Taylor series 

development of the function 



d H- 



Por small 
they wrote 


L, 


they limited themselves to the first term and 



(56) 


B]y Integration we find the "log-linear law" which Is well known 



(57) 


The problem then reduces to determining the constant a. The /25 
value 0.6 first proposed by Monln and Obukhov was found to be 
small afterwards. Among the various determinations we can mention 
those of Deacon (1962), Panofsky, Blackadar and Mac Vehll (1962), 
who found 4-4.5 and 6 for unstable conditions, respectively. 

Por stable conditions, Mac Vehll proposes the value ?• Also, 
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Zllltlnkevltch (1970) studied the various papers and gives an 
average value for a close to 1 for unstable conditions and 9 for 
stable conditions. 


The dispersion in the obtained values can be understood be- 
cause of the experimental uncertainties associated with the length 
L. In addition, the limitations in the expansion to the first 
term does not seem valid except for small values ofj •— , which 

seriously limits the range of application of the log-linear law. 


Iri' this case, the second order term comes into play and 
obviously will prevent a correct determination of a. 


This difficulty led a certain number of authors to find 
formulas which were valid over a more extended range. Among these, 
two are generally used. 


The first formula carries the name KEYPS- (formed from the 
initials of various authors which Independently proposed lt,l 
Kazanskl, Ellison, Yamamoto, Panofsky and Sellers). It can be 
derived from the turbulent energy balance using a few simplifica- 
tions. It is written i 



(58) 


and integration' leads to the following velocity profile (Klug 196?) 


with 


u^z) -- 

^ ^ r - X Arc tg r - Arc +-W.Y 
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7“' - = p - 

y'> L 

yj* . ft z. 

3 Is determined experimentally: a value close to l 8 was obtained. 


For small 


z 

LJ 


the formula of KEYPS tends toward the 


/26 


expression of Monln and Obukhov with , or 


Another formula was established by Swlnbank (1964), who 
considered a quantity X (which Is In a certain sense a generaliza- 
tion of z weighted by the thermal stratification effects) such 
that the relationship -fe. X A \ is always satisfied. The 

production of global kinetic energy Including the effect of 

Archimedes forces is written as • Swlnbank 

^z/ ' 


postulates that this quantity Is equal to 

hypothesis makes It possible to determine 
between the two relationships. He obtains 


■ax. 


This 


by eliminating X 




(59) 


which results In the following "exponential" velocity profile 

u (z) = Hjs Lo<, 

M ^ Zo 

For small j| Z ^ ^ 

We again find the formula of Monin-Obukhov with = 0 , 5 . 
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These various forms can also be discussed theoretically. The 
formula of Monln-Obukhov, If It seems valid In the range ( <^'i| , 

Is not strictly applicable outside of this Interval. The two 
other forms have a wider range of validity. However, the formula 
of Swlnbank Is also based on Intuition as well as a physical fact . 

The formula of KEYPS uses a physical law which Is better established, 
but the approximations used are rather coarse ones. Experimental 
verifications should be made to validate them. This comparison 
was carried out by Bernstein ( 1966 ) and his results are given 
In Figure 8. The formula of KEYPS seems to be best adapted as 
other authors have also stated. However, none of the formulas 
seems to have any advantage over any of the others . 


In this study, the absolute value of -| does not exceed 1. 
Under these conditions the formula of Monln-Obukhov with o<=.0,6| 
seems to be Just as valid as the formula of Swlnbank. The formula 
of KEYPS seems to be slightly closer to the experimental results, 
but the resulting velocity profile would be more complex. We 
have adapted the simplest form, that of Monln-Obukhov. 


442. For thermlal processes 


( 34 ): 




Is defined by equation /27 


k H 


( 60 ) 


According to (33) we have 




^ 4?. u » z. 


kTH 


(61) 
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We then derive: 


kH ^ 

<1)tt is therefore derived from (|)., if the ratio IT” ~ ^ H is 

n rl ■ ■'M 

known. According to equation (62) it is a function of 

kH 

The ratio ~ has been treated in many papers in addition to 
the investigations of the function (jijyj. Various and sometimes 
contradictory results have been obtained. This can be explained 
by two factors: 

1. The fact that the turbulent j fluxes t and S are simul- 
taneously and strictly conserved in the vertical direc- 
tion should be verified by experiments. This "wojul|d| 
require experimental conditions which would be hard to 
find (see Chapter II). Also there is a problem of the 
radiation flux Influence which is what is being treated 
in this paper. 

2. The turbulent] flux S must be measured directly and 
Independently. Because of the development of a correla- 
tion measurement, which has been relatively recently 
developed and is not yet well known (Buslnger etal.]l967), 
certain results were derived from the surface energy 
balance. Often they have a considerable error. 

2 

In the unstable range increases with — and can reach 

high values on the order of 3 - The upper limit is 3.5 according 

to Monin (1970). On the other hand, in the stable range, ilH 

k^ 

decreases and reaches values on the order of 0.5 - 0.6. 
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The value of the ratio under adiabatic conditions is not 
determined exactly. For a long time It was assumed that 
Recent experiments ( Zllltlnkevltch 1970 ) give values on the order 
of 1.3 which corresponds to a turbulent Prandtl number close to • 

0 . 7 . We prefer to use the classical value 1. 

As far as the exact form of the function !^li | Is concerned, /28 | 
only a few relationships have been established. For unstable 
.conditions, Swlnbank (I 969 ) proposed an empirical equation 

- which is the result of experiments, the quality 

of which was not discussed, and which corresponds well with the 

results obtained by other authors (Buslnger 1966, Record and| Cramer 

1966 , Charnock 1967 , Cyer 1967, Swlnbank and Dyer 1967, Deardorff 

1968 , Zllltlnkevltch 1970). For stable conditions, only a few 

measurements have been carried out (Record and Cramer 1966 ). A 

wind tunnel study performed by Cermak and Arya (1970) resulted In 

the relationship =C?,5‘/'3] which corresponds well with obser- 

VU/ i 

vatlons of other authors. 

Taking the continuity conditions Into account for |i| ^oL 

we therefore adc^pted the following relationships for (see 

Figure 9) : 


5 <- 0 , 02 . 


-if 

x>\ 



- 0 02. < 2 ^ < q ooi 


hh - ± 


L 


> o,ool 



(63) 


(64) 

( 65 ) 
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As far as the temperature profiles are concerned, we obtain 
the following equations by Integration of equation (48) 


Y(!f vw 

- 0,0L < ^ < 0,001 e(iL) = 6o+ 0^ Log f- 

L. Z.Q 

> o,ooi S 0(o^ooi l) 4-2 o 


^ -c - 0,02. e (2) . @ [opl L)+ 


( 66 ) 


(67) 

( 68 ) 


443 . For the water vapor transfer 



by means of the relationship 



related to 


(69) 


The problem Is Identical to the one encountered In heat /29 

transfer. The determination of the ratio -~1 Is more delicate 

K 

ku 

than determining the ratio because of the particular 

difficulties associated with the exact measurement of humidity 
and Its fluctuations (Leducq 1970). It has since been assumed 
(Dyer 1967, Swlnbank and Dyer 1967 ) that the heat and the water 
vapor are transported In a manner which Is analogous, at least 
to the first approximation. In any stability range. Consequently 
we may assume the equation 

•<H = *<E (70) 


*Translator ' s note: Equation ( 69 ) was not Included In the 

foreign text. 
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Under these conditions, the temperature and humidity profiles 
are connected by a proportionality relationship and may be 
obtained 


f(z) - f, ^ (e(z) - 6o) 
^ 2 ) 


(71) 


45 . Summary: Determination of the temperature and humidity 

profiles 

In conclusion we will specify at the beginning arbltjrary 
but plausible values of the fundamental variables u^ , L and 0^. 

The roughness parameter ' z^is derived from u* according to 
the empirical relationship of Wu (I 969 ) (Figure 7). 

If Ujf, L and 0^ are known. It is possible to calculate 
and consequently s = - . Then © (z)j can be deter- 

mined by formulas (66), ( 67 ) and (68). 

The Bowen ratio Is then calculated from the temperature 
difference using the empirical relationship reported by| 

Roll ( 1965 ). The humidity profile Is then obtained from (71) 
and J is given byS^...^. 

~i~ 

Ac J 

Thus from the three parameters u^, L and 0^, the theoretical 
profiles of temperature and humidity are directly determined as 
well as the quantities z^, 0*, B, F#, sj and J. 
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V. METHOD OF CALCULATION OF THE DIVERGENCE OF THE LONG 
WAVELENGTH RADIATION FLUX IN THE LOWER LAYERS OF THE 
ATMOSPHERE 


51 . Mechanism of the radiation flux divergence. Method of 
calculation 


If there are no liquid or solid particles (fogs, clouds, etc.) 
the divergence of the longwave radiation flux Is essentially 
caused by the presence of water vapor and carbon dioxide. These 
two gases have a spectrum which Is characterized by considerable 
absorption In the infrared range. Depending on the H 2 O and CO 2 
content, each atmospheric layer absorbs a certain fraction of 
the radiation which It receives, and radiates a certain energy, 
depending on its temperature. For thick layers, this mechanism 
can be assimilated Into a conduction phenomenon] according to 
certain authors (Brunt 1939)* On a scale considered here (between 
the surface and 10 m) , such a simplification cannot be used even 
as a first approximation. 

The evaluation of the radiation flux is primarily carried 
out on the meteorological scale. Relatively simple methods which 
allow routine determinations can be used: these are the diagrams 

of Elsasser, Kew, Yamamoto, etc.. These methods are not adopted 
for our study ,] because they apply to layers which have a thickness 
of several hundred meters and are not satisfactory on the scale 
of the surface layer. It Is necessary to use the basic radiation 
transfer equations, which were presented in a convenient way 
for this range by Brulnenberg (referenced by Kondratiev 1965 ) • 

This work was extended by Brooks (1950) and by Punk (196I) who 
developed numerical methods for estimating the divergence of the 
radiation flux In the Immediate vicinity of the ground. The 
method of Funk, more than the method of Brooks, has already been 
used by several authors (Elliott 1964, Hamilton, I 965 , Lieske and 
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Troschein 196?) and seems to be the more satisfactory one for 
the problems considered here. In order to correctly apply it, 
the mechanism which It represents must be well known. A certain 
number of precautions must be taken during the calculation. We 
analysed this In great detail. The calculations necessary to 
obtain results are long and complicated, andj numerical calculation 
must be performed with a computer. This led Gaeuskaya et al.| 

(1963) to reject the methods of Brooks and Punk because they were 
too complex. It seems that the phenomenon cannot be represented 
adequately by simple calculations. 

52 . Method of Funk (I 96 I) /3I 


We will follow the development of Punk as he presented it. 
Let w be the optical thickness measured from the bottom (see 
paragraph 53). o, r and s are subscripts which correspond to 
values] of the variables on the surface at the reference level 
and at the top of the absorbing atmosphere, respectively. 


Zs 

Zo 



Ws = O 
w 

W+dw 

Wr 

I ' 


/ ■ / - 7 7 - / / / / / / y / / / ’^0 


The radiation flux directed downwards of the layer having 
optical thickness w^ extends from the top of the atmosphere down 
to the level which Is being studied and which Is received from 
this level and Is given by (Moller 1957)j 



~o W 


( 72 ) 
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where ( 6 )| Is the emission intensity for the black body (Planck 
function) at the frequency v, corresponding to the temperature 0. 

TTy(wr - w)| is the spectral transmission of the layer between 
w^ and w at this frequency. 

is a function not only of the optical thickness but also 
of the temperature and pressure. H^O is the principal] absorber 
here (Brooks 1950, Punk 1961). The optical thickness of the water 
vapor will be substituted for the global optical thickness. The 
possible effects of CO 2 could be approximately taken Into account 
with certain corrections. 

For an Isothermal layer, the emissivlty e Is defined as the 
ratio of the total radiation flux emitted through one of Its 
boundaries to that which would be emitted if It could be considered 
as a black body. 


According to the Klrchoff law 




e Is a function of the same quantities as * 


In the mlcrometeorologlcal range, the pressure Is practlcally /32 
constant. Also It Is possible to Ignore the variation of e with 
temperature because the deviations are so small, as Deacon (1950) 
had assumed and as Staley and Jurlca (1970) have verified, e then 
only depends on the water vapor concentration (except for possible 
correction for carbon dioxide). 


*Translator ' s note: Symbol was omitted in the foreign text. 
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then becomes 



( 74 ) 


In the same way. If the adjacent layer underneath Is treated 
as an isothermal atmosphere having infinite thickness, the 
ascending flux R at the level is given by 


OO \ / 



w 


.(75)1 


0„ is the equivalent radiation temperature at the surface, 
l.e. the temperature of a black body which emits the same flux. 
At the level z , the divergence of the radiation flux is: 


hz. L 1 ^/ 3 ^ 


— r 


(76) 


or 




6r 

v\jr 


, \A)r* A'^0 

_ / A-. s^ejw-wrl 




J 


<iw| 


Wr 




(77) 


As we will see in the following, the global emlssivlty is 
not known with great accuracy. Its second derivative is there- 
fore also relatively poorly known. This is what led Brooks (1950) 
to integrate expression (77) by parts so that the first derivative 


9Q 

appears, which is better determined, and the term -c-^i 

okV 


9£i 

which is relatively well known from experiment. 
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We then obtain 



In our particular problem, the temperature and humidity /33 

profiles were known analytically. A second Integration by parts 
results In and , which was done by Punk and which 


Is of particular Interest. This makes It possible to obtain the 
following formula which will be used as the basis of our calcula- 


tion . 



53 . Definition of the modified optical density u 
The optical density Is defined by the relationship 

( 80 ) 

It Is therefore assumed positive downwards and Is zero at 
the top of the absorbing atmosphere. 


c.iw ~ p^dz. w=- 
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This definition is difficult to use in our study in which we 
are Interested in the layers close to the surface and where the 
variations in optical thickness which occur are extremely small 
in relative value. 


w does not occur in the equations except in the differential 


expressions 

"Siz. .^w 


and 


or in the form of finite differ- 


ences (Wj, - w)J . It was advantageous to replace w by its fraction, 
which varies significantly near the surface. This quantity co, 
which can be defined by the relationship uJ - 
the "modified optical thickness". 



The global atmospheric optical density is 
is therefore equal to w + coj . 



and 


CO is therefore the complement of w. Their sum is equal to 
the global atmospheric thickness. /3^ 


The various terms which occur in the formula of Funk are then 
modified in the following way 



(81) 

( 82 ) 

( 83 ) 

( 84 ) 

( 85 ) 
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The equation of Funk Is then written as 




- <rl. 


:>uj] 


^ hz 


(A eh _ 6.'* tpe^ 1 £ (tOr - Wo) 

' ° hijdr ^ W/iA^o ^ 


^ UJ UO 

f f €^(uJr-U))3u) 

S 'wib -'{Air ^ 


(86) 


54 . Calculation of the modified optical thickness 


According to our definition 


OJ ^ 


IW 


cl 2. 


^ gj Is the specific mass of the water vapor at the level z under 
consideration. If^^jls the specific mass of humid air at the 

same level, the corresponding specific humidity F Is defined by 


and consequently 




In the region under consideration, the variations of as 
a function of z are negligible to first order and 


UJ ~ 


"Pj< 


F.d 


(87) 


therefore represents a measure of the height weighted by the 
water vapor concentration In a certain sense. 


By Integrating equation (87) in parts, we obtain 
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According to (32) Z. - p 

2)Z 


4a (t-) 


or 


tx) = p 1^21 f (z.) - /4e(%)^ 


(89) 


Depending on the range of -^j considered, ai is represented by 
different analytical expressions, corresponding to the expressions 


adapted for 

*Por thermally neutral conditions, and for - 0,02 ^ ^ 0, GDI 


(90) 


*Por over-adlabatlc conditions ^ - o,02 


*Por Inversion conditions ^ 0,001 

4 >1+0,6 r V 

=> 1;(^) - 


(91) 


(92) 


55 • Calculation of the expression 


We have 


^UJ 


IS) 


«l5sl 

QuiJ 


According to (87) and, according to (31) 


c)z 


S • ¥'Wa 
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Thus , 


*For thermally neutral conditions 


"Suj ^ 


for -0,02^ - <'0,001 

*For over-adlabatlc conditions 


— - O/fe -s- 


for 


-<- 0,02 


*For Inversion conditions 


for ^>.0,001 


^ -A-+ o>^ % 


(94) 


(95) 


(96) 


56. Calculation of the expression -^rzjz 


/36 


We have 


0yr 

f 






= i. (^\ ^ ^ r_iiL 

4h 


y 




' \ _A '^<py\ 

Zf(z) “4)z 


B: 






_ 

Bz 


+- ZlfiJz- 

<^(z) ^ 




From (32) because (^h 


U„+5ti‘- zi§ 


T|hus - 
expressions 


C 


<^4^h 


which leads to the following 



*Thermally neutral conditions 


for - 0,02^ - ^D,001 




' ' ^ ml 


( 98 ) 


*Por over-adlabatlc conditions 


^ t: _ L 2S^[ 
Boy- [2.7 f ‘ 




(99) 


*Por inversion conditions 


(p2f|4 


( 100 ) 


57 • Simplifications of the modified Funk equation (equation 

86 ) 


Funk ( 1961 ) utilized a development limited to first order for 

the quantities and , taking into account the small 

Bwj 

temperature deviations which are observed in the surface layer. 


We will now do the same, utilizing the expressions which were 
established in the preceding sectionsi (55 and 56 ) in order to 
verify tliat the second order terms are Indeed negligible.' 

In effect, y ^ ( /(. 3 . ..I 


The relative variation of 0 is small (on the order of 1%) 


and we may write 


kBo — 


( 101 ) 
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According to the expressions which were established for ^ 

and for adiabatic conditions, and assuming that is small 

compared with P, which we will justify in the following, we may 
write uj (k pzf(z:)| (see (90)) and we have: 


3 ulj 


uo 


^CJU^ 




Thus 


M 

'Bwj''' 


'h. © 




and 


■d-f ^ 0: 4 'L!^ 

^ ^ ~ e 

C>IV*- 


Thus 3) 

& 


always is less than 1% and we can therefore write: 




ru 






2>u> ‘ 


(102) 


The expansion carried out by Punk therefore seems justified 
and the Punk formula can be written in the following simplified form:| 

. ■PrvrsTTt • 

=^cTBoM 2^1 |(eo-ea)^('^'-:^-gs. s£(^.-.-w) £(^,-uio) 


C>Z. 
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07 ./ 
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^(vAJs-tOf) 4- r 6l(lvu'UJrl') 1 

' UWuis JoJo J 


(.XU:}; 




|^«£-(uJr -u;o) - t [ujs-(jcr) 


'(Jc, 


^ £(uOr -UA,) cioJ ^ 

5 ! 



Formula (103) then finally becomes: 




iz.ji 






(6o-0/\') (u)r - UJj) Oq ’^(vjJS-UJt) /ci>6\ 




f U/S 


^0 

2)« 


— I S.(uJr-UJt>) - t [jiiU- UV\j( C^U> 


Wi, 


(104) 


58 . Discussion of the emlsslvlty curve (zJ}^\ 

By definition, e Is the global emlsslvlty obtained by 
Integrating over the entire frequency spectrum corresponding to 
the optical thickness (o . It can be determined either by Integrating 
the elementary emlsslvlty.- , corresponding to the various wavelengths 
of the spectrum of directly from absorption measurements of Infra- 
red radiation for various optical thicknesses. 

In the surface layer, the variation of e Is a function' of the 
pressure and temperature can be ignored (Deacon 1950) and we will 
utilize the values obtained for the atmospheric pressure on the 
ground and for a temperature close to 20° C. 

As a comparison, we studied the values of e given In the 
articles of Brooks (1950), Deacon (1950), Punk (I 96 I), Huhn (1964), 
Zdunkowskl and Johnson ( 1965 ), Atwater ( 1966 ), Stone and Manabe 
( 1968 ), Staley and Jurlca (1970) and the work of Fleagle and 
Buslnger ( 1963 ), Kondratiev (1965), Haltlner and Martin ( 1967 ). 

This revue of the literature made it possible to establish the 
following table where the emlsslvlty values obtained by various 
authors are shown for the corresponding optical thicknesses (the 
emlsslvlty values are given In %) 
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*1 


1 

1 

10-^ 

10“V 

-3 

10 

10"^ 

10 -’ 

1 

1 

1 

BROOKS (1950) 

0,29 

2,44 

11,55 

24,56 

39,59 

57,5 

DEACON (1950) 


2,8 

11,8 

25 

48,6 


HOWARD et al. (1955) 

• 

2,6 

8,9 

19,8 

36,5 

52,5 

cited by STONE and MANABE 

(1968)| 






i ELSA55ER (1 960) 


5 

14 

28 

41 

54 - 

! MOLLER and|ZDUNK0W5KI (1962) 

0,27 

2,47 

11,5 

28,9 

45,4 

59,5 

KUHN (1964) 

0,58 

4,0 

12,8 

24,9 

38,1 

54,3 

STALEY andj JURICA (1970) 

1,96 

5,65 

14,3 

27,3 

41,8 

60,0 

*,Commas represent decimal 

points .J 






Examination of Figure 10 and 

study 

of this 

table 

show that 


there 1^1 quite a large dispersion 

which 

Is explained by 

the 



difficulty of experimentally determining the function . 

Very few measurements have been made for 0)^10-^ cm H 2 D,jwhlch is 
particularly unfavorable for a study where small optical thick- 
nesses must be taken Into consideration. 

The consequences of these two remarks will be examined in 
the following. 

It seems necessary to select a more valid emlsslvlty curve. 

We agree with Kondratiev ( 1965 ) who, according to our comparative /39 
study adapted the curve proposed by Brooks D. L. (1950) (derived 
from experiments of Brooks, F. A. and Robinson). We agree with 
this because of the quality of the corresponding measurements 
and the relative Intermediate position of these values. Accord- 
ing to recent works of Kuhn (1964) and Staley and Jurlca (1970), 
it seems that greater emlsslvitles must be considered In the 
range 0 / 4 . 10 -^ cm H^O'j . It seems therefore preferable to confirm 
this tendency by additional work. 
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We would like to note that the emlsslvitles considered do 
not correspond to absorption by water vapor and therefore do not 
consider carbon dioxide. This question Is quite controversial. 
Brooks ( 1950 ) justified the omission of CO 2 by again using the Idea 
of ElsaSser according to which the quaslcomplete absorption within 
the r ■ band at the temperature under consideration leads to a 
negligible role for CO 2 In the divergence of the radiation flux. 

This point was contested by Cooley (referenced by Haltlner and 
Martin (1967), page 256.) According to Deacon (1950), the statement 
of Elsasser would be valid for high layers of the atmosphere, 
but absorption by layers at heights less than 100 m would not be 
sufficient to make the role of CO 2 effectively negligible. It 
seems more correct to take Into account the Influence of CO 2 by con- 
sidering an additional absorption which can be directly added to 
the water vapor absorption, as Deacon did. This type of calculation 
assumes that the absorption bands do not overlap except over a 
very small wavelength Interval, which seems to be valid for the 
lower layers of the atmosphere because of the small yalues| of 
optical thickness (Deacon 1950). Deacon therefore established 
emlsslvlty curves based on values of Brooks for water vapor but 
corrected by an additional emlsslvlty due to CO 2 , which Is 
assumed to be present for an average concentration of 0 . 03 ^ In 
volume. (This hypothesis Is obviously an approximate one but Is 
acceptable because It only Influences a correction term) . The 
values obtained by Deacon were utilized by Punk (196I In the 
application of his mmierlcal method. We would like to state 
that within the framework of our study, the possible effect of 
CO 2 Is of no great practical Importance, because of the disper- 
sion of the emlsslvlty values corresponding to absorption by water 
vapor only. The values of Deacon, which take Into account the 
CO 2 effect, even though they are obviously larger than those 
of Brooks, are within the fluctuation range which Is shown In 
Figure 10. The two problems cannot be separated. The only 
Important point to know Is whether the Indeterminacy of the 
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emlsslvity curve afTfects the results in a significant way. The 
calculations will then be carried out with two different emissivity 
curves 5 the one of Brooks (1950) and the one of Deacon (1950). 

By comparing the obtained results we can evaluate the role of CO2 
and the Influence of the different emissivity v|alues for water 
vapor alone at the same time. 

Another unfavorable point must be considered. This is the 
scarcity or absence of values of e for cn . In our 

study, optical thicknesses as small as 10"® cm H^O must be taken 
into consideration. Just like Zdunkowskl and Johnson (19^5) 
and Atwater (1966), an extrapolation relationship must be used. 
Assuming that e must be essentially proportional to o) (Fleagle 
and Buslnger 1963) for very small values of optical thickness, 
it seems that an equation having the form ^ (1 + o v J] 

is more approaprlate for cm H^cj (a and b are determined 

for each case t>y continuity conditions at the point 10-^ cm H 2 C| 
of |u;)| and ) . 

VI. NUMERICAL CALCULATION OF THE RADIATION FLUX DIVERGENCE: /40 
METHOD AND RESULTS 

61. Conditions for carrying out the numerical calculation 

If' we know the analytical expressions for the temperature 
profile, humidity profile, modified optical thickness m, and 
functions .8^1 and , the calculation of the radiation flux 

“3ujj 

divergence according to formula (104) can be done numerically 
for different values of u*, L and 0^, if 0^ and are defined 
and if £(u))l is known. 

Continuing the discussion in paragraph (58) and after per- 
forming numerical studies on the computer, the following expres- 
sions were obtained for representing the emissivity curves. 
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(105) 


According to 

'brooks U)^10“'^ cm H 2 O 

UJ<10”^ cm H 2 O 

According to 

DEACON U/^10“^ cm H 2 O 
j UJ^10 ^ cm H 2 O 

z was defined as the altitude at the top of the absorbing 
s 

atmosphere. This Introduces a problem, because we assume certain 
variation laws for the temperature and humidity which are obviously 
only valid in the surface layer. By continuing them and extra- 
polating them to higher altitudes, unreasonablle results are 
obtained. As Punk ( 196 I) showed, and as we were able to deter- 
mine in our calculations, the Influence of the upper layers (above 
100 m) is very small. The temperature and humidity gradients 
are primarily of Importance in the very low layers. It is there- 
fore possible to limit oneself to a value of z^ on the order of 
100 m, which amounts to assuming that the temperature and humidity 
gradients are negligible in the upper layers. 


£= 0,5777 + 0,0847 Log Ul) + 0,0027(LogUJ) 


£= 0,260 Log (1 + 'l ,04.10^ LO ) 


( 106 ) 


(107) 


£, = 0,7875 + 0,1420 Log UJ + 0,0065(Log (XJ) 


t = 0,055 Log (1 + 7,00.10^ a) ) 


( 108 ) 


Also it is possible to Ignore the terms and 

in equation (104) as Punk ( 196 I) showed. 
The validity of this approximation was verified numerically and 
can be Justified directly. In effect, the first term is close 


to 




, where o)„ is small compared with m , 
n s 


as: 

S>uj 


a rapidly decreasing function as we shall see later on. 
quently is small. Also the second term is close to 


IS 

Conse- 



Thus, £^uJ 4 ^[ is always smaller than 1 and 
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Is on the order of 7 ^j ^ which Is very small compared with the 

other terms. is large with respect to the values of w within 

the surface layer. 

The problem posed by the definition of 0^ is a more delicate 7^1 
one. intervenes through the term ” ^rI > which represents 

the deviation between the real temperature and the equivalent 
radiation temperature at the surface. 0^^ is an experimental 
quantity which is difficult to introduce into the calculation. 

Also, in the case of the ocean, whose surface emlssivlty is very 
close to 1, it seems that 0 o “ ^ r' does not exceed 0.5° C 
(Lecomte and Deschamps 1970). ' Consequently, this term is much less] 
important than above the ground where the deviation can reach 
10° C. In order to avoid the introduction of additional variables 
for the numerical calculation, we will assume that “ ^oj and 

we will later on examine the consequences of this hypothesis. 


Under these conditions, equation (104) finally becomes 



(109) 


This relatively simple expression has advantages, because 
of the form of the integral in the brackets. In effect, if we 
again consider the initial expression (79) used by Funk (1961) 
Hamilton ( 1965 ) j Lleske and Troscheln ( 1967 ) j we can see that 
the terms -ujo]|and are essentially 

of the same order of magnitude and are also much more important 
than the other terms. The result therefore is essentially given 
by the difference of these two terms. Separate estimation of 
them can lead to a considerable error, as Lleske and Troscheln 
( 1966 ) remarked. The establishment of an analytical expression 
for this difference makes it possible to limit the errors 
resulting from a numerical calculation in a significant way. 
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This Is Important to us and as far as we know. It has not yet been 
done when the Punk formula was used. 


The calculation of the integral (109) is too complex to be 
carried out analytically In the general case. Therefore we must 
carry out a numerical calculation. The calculations were first 
done by hand for particular cases. But the number of possible 
combinations of the data and the requirement for decomposing the 
integral Into a sufficient number of Intervals rapidly led us to 
using a computer. The numerical calculation was performed by 
M. R. Tomassone from the Biometry Department of the INRA (CNRZ - 
78 -Jouy-en-Josas) on an IBM 36 O- 5 O. The program which was 
established In collaboration with him Is given In the Appendix. 

It was built In order to determine the Influence of the emlssivlty 
curve on the results by comparing the values obtained with the 
curves of Deacon and Brooks. It is also used to determine the 
effect of a simplification of the profiles when the logarithmic 
relationship Is substituted for the complete expressions. 


The calculation was carried for 50 different situations 
corresponding to five values of friction velocity u^^ 


u* *» 0,08 m/a ' 
u* a 0,16 m/a 
u* a C,2B m/a 
u* a 0,5 m/a 
u„ a 0,65 m/a 


t^io 


0 m/s) 

■=%> z 

0 


10-5 

m 

(iJIO 

5 

m/a) 


<V 

5.1 0~5 

(•^10 

Ci 0 

m/s) 



0 

1 

m 

<‘^10 

£2^12 

m/a) 

•=^^0 

— 

5.10-^ 

<“10 

C^15 

m/s ) 



10-3 

m 


742 


and ten values of the Monln-Obukhov length. 

L (m) - 5, - 16, - 100, - 500, - 1000, + 1000, + 500, +100, + 16, +| 

and r( 2 ) - R^l are obtained at the levels 


The values of 
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1 mm, 5 mm, 1 cm, 5 cm, 10 cm, 50 cm, 1m, 5m, 


and ^1 0 m 


For each situation, the calculations were made using the 
following hypotheses: 


Hypothesis 1 - emlsslvlty of Deacon (for a concentration of 0 . 03 ^ 


of CO2) 

Hypothesis 2 - emlsslvlty of Brooks 

Hypothesis 3 - emlsslvlty of Brooks and logarithmic profile 
expression 


62 . Presentation and discussion of results obtained with 
Hypothesis 1 


Hypothesis 1 corresponds to the general solution of the 
problem^ because the complete expression for the profile is used 
and the role of CO2 Is taken Into account. 


41 


The results obtained for 
are given In the tables of Appendix II 
leadls to the following conclusions. 


under different situations 
Examination of them 


621. The form of the profiles always remains essentially 

similar to Itself. Figures 11 and .12 show typical profiles 
corresponding to the stable and unstable ranges. 


622 . The profile shape is very close to the one used In /43 

previous studies, and Is similar to both theoretical ones 
(Zdunkowskl and Johnson (1965), Atwater ( 1966 ), Lleske and Troscheln 
(1967) as well as experimental ones (Yamamoto and Kondo ( 1959 ), 
Hamilton (1965) and primarily those of Hlnzpeter and Heinrich 
(1969) In the range 0-10 m.) 
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The divergence In the radiation flux, no mater what Its slgn,| 
decreases rapidly with altitude and always maintains the same sign. 
However, the maximum does not correspond to the surface Itself 
but Is slightly shifted (about 1 mm to 10 cm, depending on the case). 
For the air layer Immediately adjacent to the surface (on the order 
of 1 mm), there Is an Inversion of the phenomena which correspond 
to heating when there Is cooling above It and vice versa. These 
results already obtained by Zdunkowskl and Johnson (1965) In a 
similar study and also by Godson (I 965 ) by a different method 
make It possible to explain the fact that the minimal nocturnal 
temperature often occurs not at the surface but slightly above 
It (1 mm - 10 cm, depending on the author). Numerous experimental j 
observations of this phenomena have been collected by Oke (1970), 
who himself demonstrated this fact. 

623 . Variation of the ratio ^ function of z Is 

shown In Figure 12 for various stability ranges. The Influence 
of the radiation transfer appears as an Increasing function with 
altitude. The divergence has the same sign. Thus, as Robinson 
( 1950 ) predicted and as Figure 11 shows. If the term 

becomes small above 1 m In height, it Is not possible to conclude 
that the Influence of radiation transfer Is negligible above this 
height and no statements regarding the conservation of turbulent] 
flux can be made. This conclusion has often been reached. 

624 . As we Indicated previously, the relative variation of 

this latter flux Is equal to the ratio . In the following 

table we will give the calculated values at the level z = 10 meters.* 


*For large wind velocities, Monin-Obukhov lengths such as 

L = 16. m| or 5 m lead to deviations ® 10 “ ^o| which are much 

greater than those produced In natural mechanisms and can there- 
fore not be considered plausible ones. Consequently, we did not 
give the corresponding results. 
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The relative variation of the turbulent flux resulting from 
the calculation Is sometimes very large and most often greater 
than 10 or 20^ ^ which could be established as an arbitrary limit 
for assuming that the hypothesis of flux conservation between 
0 and 10 m is valid. Nevertheless we should stress that these 
numbers do not necessarily correspond to the effective variation 
of S but Is the result of our method: If S Ts assumed to be| 

Initially constant, the variations In S given In the table can be 


6l 




















derived from it. The effective variation of S will be smaller in 
general, but must be taken into account completely in most cases.* 

625 . According to the above table, it seems that the deter- /^ 3 _ 
mining factor (for a given surface temperature and therefore given 
surface humidity) is the wind velocity. The smaller it isjjthe 
greater this effect will be. This is logical since the value of 
the turbulent flux relative to the value of the radiation flux is 
then reduced. For the same reason, the relative variation of the 
turbulent flux is greater at a given velocity in an inversion 
situation than it is in an over-adiabatic situation. This con- 
firms the generally held opinion according to which the influence 
of the radiation transfer is primarily noticeable during the 
night (and even more if the wind is weak). However, even in an 
over-adlabat ic situation or, in general, in daytime conditions, 
the effect can be noticed for weak winds, which was expressed in 
the work by Lumley and Panofsky (1964) and Munn (I966). 


63. Influence of the emlssivity 

The results described in the preceding chapter were obtained 
for emlssivity values given by Deacon. 


We can ask whether a priori . the important role played by 
small optical thicknesses in our study would not lead to consider- 
able deviations in our result. In effect, if the values of e 
adopted by Brooks and Deacon, for example, only differed by 10 to 
20 % for 'UD > lO""* cm H2 Oj , a much larger difference would be 
introduced in the extrapolajtlon formulas. The ratio 

DlACON 

would go to 0.67 when u goes to zero. 


*Thls point will be discussed in the conclusion. 
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The study of the ratios BROOKS and j^ %)' ^BROOKS 

r^^^D^AcoU CBI^tMpeacon 

especially for z « 10 mj j makes it possible to study the deviation 
which could result by adopting any kind of emlssivlty curve. 

These ratios, calculated for a certain number of cases, are given 
in the following tables. 



‘^o^brooks 

)deacon 



The differences in the divergence are on the order of 20^. 
The differences in the flux variation are about 10%. These num- 


bers can be considered to show the Influence of CO 2 or the un- 
certainty connected with the emisslvlty value for water vapor 


only. It therefore seems that the high values of £Hl 


obtained in the preceding paragraph cannot be attributed to 
excessive values of e,' because the contribution of CO 2 does not 
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exceed 10 %. On the other hand, the Increased einlssivlty values 
which were the result of recent studies by Kuhn and by Staley and 
Jurlca could not lead to a higher number for the relative variation 
in turbulent flux. 


In conclusion^ the problem of exactly defining the emisslvlty 
values does not seem to be the determining one, as Zdunkowskl and 
Johnson (1965) and Jurlca (I 966 ) found. Thus the latter concluded}] 
"T|hus, we may consider that the Brooks method shows that there 
Is a moderate sensitivity to emisslvlty values when considered as 
external parameters. It seems that any calculation based on a 
reasonable collection of values, no matter from what source, will 
lead to results which can be considered equivalent for all practi- 
cal applications." 

64 . Influence of stability. Logarithmic approximation of the 
profiles . 


The Influence of stability was demonstrated by the numerical 

results given In paragraph 62. The appearance of the curves for 
R(z)" ^0.1 

5 _ and those shown In Figure (12) suggests the definition 


of an empirical stability function ^ 


such that 


- .Jjp - hr<rel 0^ 

5o I 


The dispersion of the various values of 
relatively small, which justifies the method. The 
values are given In the following table 



obtained Is 
averaged 



( ^ 
I L 

. 1 - 0,6 - 0,3 - 0,1 - 0,05 - 0,02 0 + 0,02 + 0,05 + 0,1 + 0,3 + 0 , 6+1 

. 5 

1 

],89 0,91 0,95 0,99 1 1 1 1,05 1 ,06 1,10 1 ,12 1 ,25 1 , 4 , 


*Commas represent decimal points.] 
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L = 10 m is a more modern example and shows that the Influence 
of stability in the first 10 meters of the atmosphere will be less 
than 155? in the over-adiabatic regime and H0% in the inversion 
regime . 

It seems that the simplified logarithmic form, of the profiles 
will not bring about Important errors in estimating the relative 
variation of turbulent flux. The second order effects can be taken 
into account by the following empirical relationships: 


V(f) f’l for L < o| 

(110) 

+ 0,4 i)j for L >0j 

(111) 


In order to confirm this opinion we carried out the calculation 
using hypothesis 3, the Brooks emlsslvlty and the logarithmic 

expression for the profiles. Comparison with the results 

* 

obtained for hypothesis 2 lead to the following table.] f 


_z 

L 


- 1 - 0,6 - 0,3 -0,1 -0,05 -0,02 0 
3,90 0,93 0,96 0,98 1,02 1 1 


*Commas represent decimal points, 


+0,02 +0,05 + 0,1 +0,3 +0,6 1 

1 1,03 1,04 1,10 1,27 1,52 


Agreement with the preceding table is very satisfactory. 


Nevertheless, it is obvious that the loss of accuracy 
corresponding to the simplified logarithmic relationships in the 
numerical calculation is not an advantage. The reduction in 
calculation time on the computer is not significant. On the other 
hand, it becomes difficult to find an analytical solution because 
of the complexity of the relationships used. However it is 
possible if we limit ourselves to logarithmic expressions, as we 
will see. 
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VII. ANALYTICAL CALCULATION AND ESTABLISHMENT OF SIMPLIFIED /48 


RELATIONSHIPS 


71 . Interest In an analytical calculation 

Even though the calculation cannot be done except by using 
a certain number of approximations, it is of interest for two 
reasons;| 

— on the one hand, it makes it possible to have approximate 
calculation methods for the divergencej which are simple 
and fast (the computer is replaced by the slide rule) which 
are then of Interest for the experimenter 

— on the other hand, the relative influence of mlcrometeoro- 
loglcal parameters can be shown. 

7 2 . Analytical expression of- the functions 0 and m 

The temperature and humidity profiles are reduced to the 
logarithmic form 


©(z) -e. = e^LooB 

f^(z) - ~ 

o 


( 112 ) 

(113) 


According to paragraphs 54, 55, and 56, 
written as : 


U) / 


^ pj- 


are 



( 114 ) 
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( 115 ) 


f(4 


(116) 


According to the calculated humidity profiles. 


F* 

F(?) 


always /49 


remains less than 0.1. The relatively large values of this ratio 
correspond to relatively great heights. We can therefore Ignore 


the term 


ir*‘ 

fTT) 

following equation: 


In equations ll4, 115, II 6 which leads to the 



(117) 


- Q »r 
la) 


(118) 


U)^ 


(119) 


73 • Analytical representation of the emlsslvlty curve 

Such a representation must approximate the experimental 
curve In a satisfactory way at least over the optical thickness 
range used In the calculation. Also it must have a simple form 
which will not cause too much complexity In the analytic calcula- 
tion of the Integral In equation (104). 
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Such expressions have already been used for numerical calcula- 
tions. Thus, Hamilton (1965) proposes a rather simple decomposition 
of the curve £.{u)),|lnto two parts of straight lines, which meet 
at the value U) = 8.1 cm H2Q.I Atwater (I966) utilized the de- 
composition of Elliott and Stevens In which the curve Is repre- 
sented by line segments for the Intervals limited by the values 
10 “^, 10 “^ sj 10“^1 and 10“2 cm Finally, Jurlca (1966) also 

utilized analytical expressions for which enter the formula 

of Brooks,|but he does not give them explicitly. 


Even though the expression utilized by Atwater seems to have 
the required accuracy, the decomposition Into several Intervals 
Is cumbersome for the analytical calculation. We also look for 
a relationship better suited for our study. Since we could not 
find a simpler relationship which extends over the entire variation 
range of w, we had to consider two Intervals. 

Ppr lU.^10”^ cm H^O I 


the form of the emlsslvlty curve of Brooks In the seml-logarlthmlc 
representation led to the use of a hyperbolic tangent function of 
the form: 


IB 

a. IS.)" ^ 


( 120 ) 


e^, A and B are numerical constants. 

Such a simple expression cannot cover the entire variation /50 
range of o) with the desired accuracy. The Introduction of higher 
order terms would lead to too much complexity. 
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The contribution of layers such as cm H 2 O 
quite small. ¥e determined A and B such that exp] 
approximates the experimental curve over the Interval 


L cm H 2 O ^ ^°”^cm H 2 ° 


1-I 


In the best way. 


The adopted values are 

( A = 1,125 
B = 0,25 
' ^ = 0,210 


As Figure 13 shows, agreement with the experimental data is 
satisfactory over the interval Z, 10“^ cm H^O, 4 x 10“^ cm H20__y |.i 

There Is a slight divergence which appears above 4 x 10"2 cm H20I 

This Influence Is less than \%. All attempts to make an Improve- 
ment on the side of large optical thicknesses led to a divergence 
on the small optical thickness side, which was much more detri- 
mental to the calculation accuracy. 


The corresponding analytical representation of the first 
derivative Is shown in Figure l4. It also seems to be 

satisfactory considering the dispersion of the experimental 


values (the high values of 


5 >€ . 

■£>uJ 


proposed by Fleagle and Buslnger 


(1963) and which were contested by Staley and Jurlca (1970) were 
not Included). 


For uO^IO"^ cm H a I used an extrapolation formula of the 
form ^ = a Log (1 + buJ ) (Cf ^8),;|a and b are determined from 
continuity at the point 10-4 cm H 0.1 


We have 


£ = 0,260 Log (1+1,04 x 1 0^ OJ ) 


( 121 ) 
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O/C./ 

This expression as well as the one corresponding to ^ are 




also shown In Figures 13 and l4 


74. Details of the calculation proper 


751 


This calculation Is based on the equation (104) In which the 
terms ^ and Qe:(uJ. can be neglected 
as we saw. We therefore find: 

V O^m / 


Zr 
/C0$ 






( 122 ) 


Taking 117, H 8 and 119 into account we obtain: 


al.' 


( 01 = - ^“^ 7 ^ 41 '“-^) 

^ j J 


1 




Vo 






( 123 ) 


The calculation amounts to calculating the expression 

U/9 ' 


X = 






which Is done In different ways depending 


on the value of 


731 . In the most complex case oJr^1,1 x 10"^ cm H 2 O, 
must be considered. 


5 Intervals 


I. o* 


Let US: set 


jcOr -wj = 


U)r-UO 
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and UJr > u)( 'w::^ 0,6 uJr 


This latter condition leads to assuming that 


Is approximately constant and equal to • 


<S>S I UJr - U >) 
3uj 


By Integrating by parts X -1 = X (w<j ^ uji) = I 4 ^ 

_ _ a£(u)r ) n oMy^ ■ JuJo 

Slo L -3 UJo^ I 


10 


-4 


'J 

1 u)r -IXI \ ~ u)r - UJ 

eq 1 ( 120 )| 


)<i = X (o , olr - =^2 ft c, ■ (S. \1S,) 

^ ^io->+CT 3 j- ^A.al,_.A,al,^(^,j 

l^Wr - ^ 

\ 

jlkJr- u;] = UJr - td 
ujf - tu s io ^ ^ 


7 /52 


X-. = X («,-!.• , u;r) = Uj 


UJr 


Wr-'O'V 


IV. 


l^lOv' ; LxJr + -^0 

U)r - l>o| = 0 J-U)r 


UOr - UJ $ \0 


-V 




X^=x.(«.,w.+^»-‘^) 7 -/^] 

/uJr 
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V. 



1 


, UJsJ 


We are then referred to case II with j^j. 


oJ 




and we 


find : 


-1 


( 




- - 2 A£\ V UU, ^ j — I , / 

\ \ l0o-'.0-r A )JZh _ Arri?i 1 




Thus for the numerical calculation, we see that Is a 
quantity which Is not well defined. 

For the analytical calculation, we may set LUg = 100 vcyj and 
Ignore the contribution of the upper layers. 


We finally find: 



Loq '^O - Loq -iii!!. + 



g)y 
L'/r + li*- 



-s- Xs 


But and are 
depend on co^. 


defined numerical functions which always 


/53 


We can therefore set 


X 





(124) 


R^ju;r)| Is a function which can be e'xpresse'd In' analyt'lcal form 
and which can be calculated once and for all. 
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732 . In the case 10"^^ UJr ^ 1,1' . 10“^j 

0,1 UJt^uJr- 

The Intervals II and III are regrouped Into the single inter- 
val \_o,1 uJr ^ • The Intervals I, IV and V remain unchanged. 

We then obtain an expression having the form (124 ),1 Is 

this time defined by a different analytical expression. 

733 . In the case u)g^ ui r ^ Intervals I^ 

II, and III are regrouped into the interval cur'll and the 
Intervals IV and V remain unchanged. Using the same calculation 
method, we again find expression (124). and even here we find an 
analytical expression for ,R|(uJr)j • 

734 . In conclusion, in the general case and for any oj we 
can write 



Ri[u)i^j is a defined function and its analytical expression 
differs, [depending upon the interval under consideration. But 
obviously it remains continuous over the entire variation Interval 
of co^. It was numerically calculated and is shown in Figure 15 
(it should be noted that is a negative quantity for 

cm H20)| . 


735 • According to equation (123) and the expression for X, /54 
we finally find: 
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As soon as ^ 10j j which occurs most frequently, 

and we find: 

. I 



With this expression we can directly calculate from 

values of and . This Is done using Figures 15 and 

l6. The calculation Is slmple,|because It Is sufficient to 
calculate w □ = ^ z^j and u)r = P '^(zr) ^r.l 


It would also be necessary to know , which we 

Ignored In the numerical calculation. We can determine its In- 
fluence by means of an analytical calculation. Let us consider 


that : 



lOo J ^ Zo 



W'p can see that In equation (124) t),-,- 0,;^| plays the role which 

Is similar to the role of e(zr)-6o] • Its Influence Is negligible 

If this is not true and because Is still 

positive, this term tends to increase the value of for 

, l.e.,iln an Inversion situation or, on the other 
hand, to be reduced for - 0^ si Oj , l.e.,)ln an over-adlabatlc 
condition. This was confirmed by the numerical results obtained 
by Zdunkowski and Johnson ( 1965 ) and Zdunkowskl, Henderson and Hales 
( 1966 ). In any case, the small value of this term in the ocean 
makes It possible to Ignore it, except If the temperature difference 
between the air and the water Is small. In this case it would 
have to be determined by experiment. The correct determination of 
the surface temperature of the ocean 0^ can never be done except 
by means of a radiometric method, which will then result in a value 
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corresponding to 0j^. This Is lnconvenlentj|because the sum of the 
terms In parentheses In (124) Is equal to &(gr)- ff g and therefore 

o 

0Q IS not Involved. 

736 . The analytical expression which was obtained for /55 

makes It possible to establish the expression for the radiation 
flux variation R^ 2 ) - Ro| by means of an Integration: 



-Ro 

! 

S(cur)l Is a well defined function of cO j,j . This Is also 
true for and, consequently, for the Integral 

which can be called R 2 ('*^r)| • calculated 

numerically and Is shown In Figure I 6 . ■ • 

We then finally obtain: 

( 126 ) 

/56 

The analytical calculation leads to simplified formulas (125) 
and ( 126 ) which make It possible to rapidly calculate the values 




7 5l- Discussion of the simplified formulas 



75 




of and R/ N - R ( at a level z without a computer, if Uj^, 0 , 

o| 

L are] known In our problem or if u*, tfo# L, z^, a- r| and B are known 
in the general case (for example,) above the ground). These formulas 
assume that the form of the various profiles can be considered 
logarithmic. We have seen that the errors introduced by a strong 
stratification can reach in the inversion range and 20 % in the 

over-adlabatlc range. 

Considering the large number of approximations which we have 
introduced, we must now determine the error of the calculation. 

The results of the analytical calculation and the numerical cal- 
culation (both carried out for hypothesis 3: Brooks emlsslvlty 

and logarithmic profiles) wjere compared for different situations. 

In all the cases studied, this comparison leads to conclusions 
which are similar to those which can be derived from the typical 
cases shown in Figures 17 and l8. According to Figure 17 which 
shows the radiation flux divergence, the agreement seems satis- 
factory overall, even though the numerical calculation Introduces 
a certain dispersion In the points, which can probably be attri- 
buted to calculation of the Integral which occurs In the expression 
for — . The ordinates are plotted on a logarithmic scale and 

the deviation between the two methods of calculation between 1 
and 10 m, which do not look very great according to Figure 17, 
leads to an appreciable difference for the Integrated term 

_ R I j as Figure l8 shows. We are now faced with a problem 
of determining which of the two methods Is more reliable. Con- 
sidering the approximations made for the two cases and the details 
of the two methods, we believe that the numerical calculation Is 
probably afflicted with the greatest amount of errors. 

This point of view confirms the difficulty In estimating the 
radiation flux divergence accurately using the numerical method 
of Funk from a limited number of temperature and humidity 
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measurements. This also explains the general disagreement observed 
(Punk 1961, Hamilton 1965, Lleske and Troscheln 1966) between the 
direct determinations of this divergence and the calculated pre- 
dictions . 

75I. Examination of the Influence of various mlcrometeorologl - j 
cal parameters 


or 


Assuming that bo = ^ r > equation ( 126 ) leads to: 


1 




(T6o‘ 

U; 




■-:>o 
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( 127 )| 


It follows that the relative variation ojf the turbulent flux 
resulting from the effect of radiation transfer: 


- Is proportional to 

- Is Inversely proportional to 


- Increases with , l.e. with average air humidity 

- decreases slightly with 

For the ocean, the variation range of 0 ^ Is quite small and 
therefore the term 0 o^j will vary slightly. Also, z^ depends on 
u*. Thus, for the ocean, |^:^|| is practically Inversely pro- 
portional to Ujf and Increases with the average air humidity. This 
confirms the results of the numerical calculation. Also, we are 
now In a position to draw another conclusion by estimating the 
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relative effects of stability. This made It possible to draw 
Figure 19 which clearly shows the Influence of the various para- 
meters. In particular, we can see that the air Increases very 
rapidly when the friction velocity drops below 10 cm/second. 

VIII. Conclusion 

A rigorous study of heat transfer In the surface layer of the 
atmosphere requires knowledge of _the] Interaction mechanisms between 
the radiation transfer and the turbulent exchanges. Up to now, 
this topic has only been studied in a limited number of works. Its 
Influence was placed in perspective In certain theoretical and 
experimental studies. The classical theories onj turbulent heat 
transfer In the surface layers do not take into account the possible 
effects of radiation transfer. Consequently errors can be made 
either during the establishment of seml-emplrlcal predictions or 
when experimental results are Interpreted. 

We attempted to determine these errors by using absurd reason- 
ing: the turbulent flux Is Initially assumed to be constant,! and 

we then calculate the variation of the same flux which Is a logical 
consequence of this hypothesis. The numerical results show that 
In certain situations , the relative variations are much larger 
than the limit of 10 to 205? which could be established, considering 
the tolerable error levels in micrometeorology. It therefore seems 
that In the general case and a priori it Is not allowable to assume 
that the Influence of radiation transfer Is negligible. 

The dominating factor Is the wind velocity. Small values 
of It lead to large variations. Thermal stratification only enters 
in on the second order level. This point allows one to find an 
approximate analytical solution to the problem using temperature 
and humidity profiles which obey the classical logarithmic re- 
lationships. Simplified formulas were established with which It 
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is possible to rapidly evaluate the effects of the radiation flux. 

In addition, this method has the advantage of clearly demonstrating 
the role of various parameters. ¥e finish with a quantitative 
determination of the unfavorable effects, such as small wind velo- 
cities, increased humidity values and nocturnal Inversion regimes. 

In addition, appreciable errors can be made if we deliberately 
ignore the influence of radiation transfer during the day and 
during periods of low wind velocity and large humidity. 

The calculations were carried out assuming clear skies. The /59 
fact that there was no fog could only Increase the Influence of 
radiation transfer (Zdunkowski et al.|(1966) concluded that there 
was a smaller fog effect than the one we predicted). A large 
difference would lead to results which would be quite 
different from the ones we obtained. The surface emlsslvlty must 
therefore be taken into account. 


Except for these points., the uncertainty, in the emlsslvlty 
values of water vapor for small optical thicknesses (less than 
.-4 


10”** cm H 20 )j and the uncertainty in the carbon dioxide effect, 
as well as deviations observed between the result obtained 
numerically and analytically, contributed to an appreciable un- 
certainty in our results, which could be reduced by additional 
work. Even though this uncertainty limits the accuracy of the 
preceding results, it is], nevertheless not large enough to 
Invalidate their significance. 


Considered from a certain point of view, the preceding work 
can be considered to have a negative quality. This is because 
we are only able to define conditions under which the radiation 
transfer can be considered negligible. We are not able to deter- 
mine the corrections which must be made for taking it into account. 
In order to do this, a method of successive approximations would 
have to be developed. We assumed|that the flux was constanti] and 
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we ended with a variation of %j . We would have to perform a 
new calculation based on this variation of %,\ which would then 
lead to a resulting variation of X2 • Assuming that the cal- 
culation converges (very close tests make It possible to believe 
that this will be true)jjwe should very rapidly arrive at a rigor- 
ous solution of the problem. Unfortunately, this method, although 
promising cannot be developed because of the present state of 
the theory. In , fact the turbulent flux S Is not constantj|and the 
Monln-Obukhov length L must be considered as a function of alti- 
tude and therefore cannot play the role of a length scale which 
would be developed over the entire surface layer. Anj extension 
of the present theories which consider the turbulent fluxes of 
momentum and water vapor as being constant but use variable turbu~ 
lent heat fluxes, would be necessary. 

Assuming that this problem Is solved', we can assume that the 
method under consideration would converge rather quickly. In 
effect, the divergence of the radiation flux Is primarily deter- 
mined by the shape of the humidity and temperature profiles very 
close to the surface. The latter would not be greatly modified 
by variations of the turbulent flux, which ajffect primarily the 
upper layers of the surface layer. Considering the uncertainty 
which affects the calculations of the radiation flux divergence, 
we can assume that the margin of error on this latter quantity 
would very quickly rise above the significant margin for 
Iteration processes. For temperature and humidity profiles, the 
difference would not be greater than the one which results for 
velocity profiles when a constant friction term Is adopted (flat 
plate case) or when a variable 'friction Is adopted (tube case). 

On the other hand, the effect of radiation transfer will certainly 
be Important during a calculation of the values of turbulent heat 
transfer coefficients from experimental profiles In the surface 

layer. This explains In part the dispersion of the results for 

k H I 

the ratio . In order to evaluate It precisely. It would be 

kMl 


necessary either to remove the Influence of radiation transfer 
(which Is possible in the wind tunnel) or to at least reduce It 
by avoiding determinations when the wind is too slow and the air 
humidity is too great. Also the fluxes would have to be measured 
as close as possible to the surface. 

This study does not pretend to entirely resolve the problem 
posed. Instead we give a first approach which should lead the 
way to further work. This work should take into account certain 
aspects which we have separated out due to the definition of our 
problem, in particular the effects of advectlon and the effects 
of non-steady conditions, which must be considered in the general 
case . 
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APPENDIX I 

PROGRAMME DU CALCUL ,MUiMERIOUE( FORTRAN IV) 


i dimension U(5) ,20M(5) ,OBUKO(10) 

dimension 2(15),X(15),TETA(15),EHUM(15),W(15),T2W2(15),V(15),Y(15) 

,, dimension WR(10),ZR(10),DIV(10),FI(10),CUMUL(10),T3W3(15),EMIS(15) 

1 F0RMAT(////,3X, ‘VITESSE DU VENT A 10 M., U10= *,F 9.1 , • . . . .RUGOSITE, 

; ^ Z0= ',E 9.3, » . .L.MONIN-OBUKHOV L= »,F9.0) 

2 FORMAT (/ ,1H , ' .USTAR= ■ ' , FI 0 .4, » , FLUX QH*’ ^F1 0.4, • ,EVAPO= 

, 1',F10.5,',B0WEN= ',F10,4,» ,‘A= ,F10.5) 

;3 FORMATdH ,100( '*' ) ) • 

.4 F0RMAT(1H , HAUTEUR * TEMPERATURE * HUMIDITE 

,1 W * D2T/DW2 *•) 

:5 FORMATdH ,5(E15.0,4X,1H*))_ 

;6 F0RMAT(1 HO, *NP-R * WR * ’ ZR * INTEGRALE * 

■ 2DIVERGENCE * CUMUL *•) 

7 FORMATdH, 14 , 1 X, , 5 (El 3 . 7, 1 X, »* » ) ) 


INITIALISATION DES VALEUR5 C0NSTANTE5 DE'l'ANALYSE 

.,i N=14 
, AK=0,4 
‘ TET0=293. 

U(1)=2.0 
U(2)=5. 

U(3)=0. 

: U(4)=12. 

' U(5)=15. 

Z0M(1 )=0. 00001 
Z0M(2)=0. 00005 
Z0M(3)=0.0001 
i Z0M(4)=0.0005 
Z0M(5)=0.001 . 

■ OBUKOd )=-1000. 

: QBUK0(2)=-500. 

OBUKO(3)=-100. 

; 0BUK0(4)=:-16 

I 0BUK0(5)=-5 
DO 100 1=6,10 
J=11-I 

OBUKOd )=-OBUKO(J) 

E0=23.27 
A1 =13000 

MX EST LE NUMERO DE L'UNITE DE SORTIE 

8,«x=3 ■ 





VALEUR STANDARD DES HAUTEURS 

X(1)=0. 

X(2)=0. 00001 
X(3)=0. 00005 
X(4)=0.0001 
X(5)=0.0005 

X(6)=0.001 ■ 

X(7)=0.005 

X(6)=0.01 ■ ■ 

X(9)=0.05 
X(10)=0.1 
X(11)=0.5 
X(12)=1. - 

X(13)=5. 

X(14)=10. 

X(15)=100. 



VALEURS DES PARAMETRES DE LA COURSE D'EMISSIVITE 

B0=0.7675 , . 

B1 =0.1 420 
B2=0.0065 



J'INDIQUE ICI LE NOMBRE DE C0UPLES(U1 0,Z0)ETUDIE5,A L’EXPLOITATION,NTEST= 
DO 130 1=1,5' 

U10=U(I) 

Z0=Z0M(I) 


ON DOIT D»ABORD CALCULER LE NOMBRE DE POINTS DU PRDFIL ET LES 
VALEURS POSSIBLES DES HAUTEURS, JAMAIS ZERO... 

NP=N+1-I 

DO 120 J=1 ,NP ' 

K=I+J 
0 Z(J)=X(K) 


ON EFFECTUE LES CALCULS POUR DES VALEURS FIXEES DE LA LONGUEUR 
DE MONIN-OBUKHOV 

DO 130 J=1 ,10 , . • ' . . . . i 

IH=1 ' • 

DO 130 U=1 ,2 ' . 

AL=0BUK0(J) ■ ■ • ‘ V 

WRITE(3,1 )U10,Z0,AL 

UET0I=USTAR(AK,U10,AL,Z0,INDIC) 

CALCUL DU PROFIL DE TEMPERATURE 


CALL protp(al,ak,uetoi,zo,teto,np,z,teta,ix,ij) 

TET10=TETA(NP) 


CALCUL DU RAPPORT DE BOWEN 



B«B0WEN(TET1O,TET0) 
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? CALCUL DU FLUX DE CHALEUR ET DE L'EVAPORATION 

> ' CALCUL DU PROFIL D»HUMIDITE 

'I CALL FLUEV(UET0I,AL,B„QH0,EVAP0) 
call PR0HU(E0,B,NP,TETA,TET0,EHUM) 

CALCUL DES ELEMENTS INTERVENANT DANS LA DIVERGENCE DU FLUX 
RADIATIF 

A2=TET0*UET0I*UET0I/(AK*AK*AL*9.81 ) 

CALCUL DES DEUX TABLEAUX CONTENANT W ET LA DERIVEE SECONDE 
DE TETA PAR RAPPORT AW. ' . ' 

call EOPTKA1 ,A2,B,NP,EHUM,Z,W,T2W2,AL,IJ) 

WRITE(MX,2)UET0I,QH0,EVAP0,B,A2 

WRITE(MX,3) . , , . . 

WRITE(MX,4) , , 

WRITE (MX, 3) 

DO 135 K=1 ,NP 

35 WRITE (MX, 5)Z(K), TETA (K),EHUM(K),W{K),T2W2(K) 

CALCUL DES VALEURS ZR ET WR 

DO 140 K*1 ,10 , * 

L=NP-K+1 . ' 

M=11-K 

WR(M)=W(L) - . ■ ' 

JO ZR(M)=Z(L) . , , , , 

CALCUL DES EMISSIVITES . 

• > 

38 DO 1 50 M=1 ,10 
, DO 1 63 K=1 ,MP 
V(K)=W(K)-WR(M) 

IF(V(K))164,162,164 
52 EMI5(K)=0. 

GO TO 163 

54 IF (AB5(V(K))-1 .E-4) 165,165,161 
51 IF(IH.NE,3) GO TO 300 
IX=2 

EMI5(K)=0,218*(1 .+1 .125*(50RT(ABS(V(K)))-7.9E-2)/(50RT(ABS(V(K)))+17.9E-2)) 

ED TO 163 

' i5 IF(IH.NE.3) GO TO 301 
IX*2 

EMI5(K)=0.260*AL0G(1 .+1 .043E3*(ABS(V(K)))) 

GO TO 163 

11 EMI5(K)=0.055*AL0G(1 .+7.E3*ABS(V(K))) 

13 T3W3(K)oT2W2(K)*(EMIS(1 )-EMIS(K)) - 

INTEGRATION DE LA '’ONCTION T3W3 TABULEE DANS L'INTERVALLE V(1) A V(NP) 
L» INTEGRALS VAUT FI.SSP PAGE 289 
84cALL QTFG(V,T3W3,Y,NP) 

FKM)-Y(NP) 



j CALCUL DE LA DIVERGENCE AU NIVEAU ZR 

I L=NP-M+1 

jfl DIV(M)=-22.4E-9*TETQ**3*EHUM(L)*FI(M)/A1 

I CALCUL DE5 AIRES CUMULEE5 DE DIV EN FUNCTION DE Z 

i; CUMUL(1 )=1 .E-3*DIV(1 ) 

DO 170 M=2,10 ' 

!o CUMUL{M)=CUMUL(M-1 ) + (ZR(M)-ZR(M-1 ) )*(DIV(M)+DIV(M-1 ) )*.5 
^ WRITE(MX,6> 

I DO 180 M=1,10 ' 

:0 WRITE(MX,7)M,WR(M),ZR(M),FI{M),DIV(M),CUMUL(M) 

I GO TO (210,130), IX . . . 

io IH=3 • . ■ • - ' 

j GO TO 208 . • . 

;C CONTINUE 

: CALL EXIT , ' . 

I END . 

SOUS-PROGRAMMES 


: FUNCTION U5TAR(AK,U10,AL,Z0, INDIO ^ 

CALCUL DE LA VITESSE DE FROTTEMENT USTAR{U ETOILE) 

. .-.AK. . .CONSTANTE DE KARMAN 

, .. .U10. .VITESSE DU VENT A 10 METRES .. ' 

.. .AL. . .LONGUEUR DE MONIN-OBUKHOV 
. . .ZO. . .PARAMETRE DE RUG05ITE 

...INDIC..0 SI LE CALCUL EST NUMERIQUEMENT POSSIBLE 

..1 SI LE TERME INTERVENANT DANS LE LOGARITHME EST, NUL ON POURRA AMELIORER 
CE POINT EN IMPOSANT AU TERME D’ETRE 5UPERIEUR A EPS 

INDIC=0 ' ' . . 'v, .■ ’ ■ 

; EPS=0. ■ ■ 

j X=AL*(EXP(10,/AL)-1.)/Z0 
: IF(X-EPS)100,100,110 

^0 INDIC=1 
■ USTAR=0. 
i RETURN 

f U5TAR=AK*U10/AL0G(X) ' 

RETURN ‘ ■ ' ' ' 

ENi ■ ' 
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i subroutine PROTP (AL,AK,U5TAR,ZQ,TET0,N,Z,TETA,IX,IJ) 

1 DIMENSION Z(1 ),TETA(1) 

I 

I CALCUL D'UN PROFIL DE TEMPERATURE TETA(Z) EN FONCTION DE Z 

tableau CGNTENANT LES DIFFERENTES VALEUR5 DE LA HAUTEUR DES OBSERVATIONS 

1 dimension N 

i .. .AL. .. .LONGUEUR DE MONIN-OBUKHOV. 

I ...AK....CONSTANTE DE KARMAN. 

j ...USTAR. VITESSE DE FROTTEMENT 

I ...ZO PARAMETRE DE RUGOSITE 

! ...N NOMBRE DE VALEURS DIFFERENTES DE LA HAUTEUR 

...TETA.. TABLEAU CONTENANT LES RE5ULTATS DE DIMENSION N 
. ..TETO..VALEUR FIXE DU PROFIL POUR 'ZO 
GO TO (201 ,202), I J 

T=ABS(AL*EXP(-0.02)-1 .) . » 

ix=i ■ ■ . ;■ 

G=9.B1 

T=TET0*(1 .+U5TAR»USTAR*ALOG(T/ZO)/AK/AK/G/AL) 

DO 90 1=1 ,N 
IF(AL)100,100,121 

lO BL=-.02*AL • ■ . 

IF(ABS(Z(I)-BL) .EL.0.000001 )Z(I)=BL 
IF{Z(I)-BL)110,11D,120 ■ 

CL=ABS(BL) 

TETA(I)=T+2.8*USTAR*USTAR*(1 ./Z(I)»*.5-1 ./CL** .5)/G/ABS(AL)*^^ .5 ■ 

GO TO 90 

S=ABS(AL*(EXP(Z(I)/AL)-1 .)) 

IF (ABS(S) .L£. 0.00001 )TETA(I)=TET0 
IF(TETA(I).EQ.TET0)G0 TO 90 

TETA(I)=TET0*(1 .+USTAR*USTAR*ALOG(S/ZO)/AK/AK/G/AL) 

GO TO 90 

DL=0.001*AL ■ . 

IF(AB5(Z(I)-.DL) .LE.0.000001 )Z(I)=DL 
IF(Z(I)-DL) 110,110,122 
TPRIM=AB5(AL*(EXP(Q.001 )-1 )) 

TPRQM=TET0*(1 .+USTAR*USTAR*ALOG(TPRIM/ZD)/AK/AK/G/AL) 

TETA(I)=TPROM220.*USTAR*USTAR*TETO*((Z(I)/AL)**0.1-O.D01**0.1)/AK/AK/G/AL 

0 CONTINUE . ^ 

O RETURN ■ ■ ' 

TETA(I)=TET0*(1 .|+USTAR**2*AL0G(Z(I)/Z0)/AK/AL/AK/G) • ■ : • 


■y-- 


IX=2 
RETURN 
END 
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FUNCTION B0WEN(TET10,TET0) 

CALCUL DU RAPPORT DE BOWEN 

DELTA=TET10-TET0 
IF(DELTA)100,100,110 
,|Q B0WEN=.1-.09*DELTA 
MX=3 
RETURN 

BOWEN=-.0*DELTA 
RETURN 
END- 



SUBROUTINE PLUEV (US TAR, AL, BOWEN, QHO.EVAPO) 


CETTE SUDROUTINE CALCULE 

..QHO FLUX DE CHALEUR SENSIBLE i ^ 

. .EVAPO. EVAPORATION b 

C0NNAIS5ANT 

. .USTAR. VITESSE DE FROTTEMENT 

..Al .LONGUEUR DE MONIN-OBUKHOV.NON NULLE 

. .BOWEN. RAPPORT DE BOWEN, NON NUL 

QH0=-.9E+04*U5TAR*USTAR*USTAR/AU • 

EVAPO=1 .5E-02*QH0/B0WEN 

MX=3 

RETURN 

END 


, SUBROUTINE PROHU ( EO , BOWEN, N, TETA , TETO,EHUM) 

; DIMENSION TETAd ) ,EHUM(1 ) 

CALCUL DU PRCFIL D'HUMIDITE C0NNAI55ANT CELUI DE TEMPERATURE 

...EO VALEUR DE L'HUMIDITE POUR Z=ZO PARAMETRE DE RUG05ITE 

i ...BOWCN.. RAPPORT DE BOWEN 

...N NOMBRE DE VALEURS DIFFERENTES DE LA HAUTEUR 

...TEtA... TABLEAU DE DIMENSION N CONTENANT LE5 VALEURS DU PROFIL DE TEMPERATURE 
...TET0...C0N5TANTE CORRESPONDANT A LA VALEUR DE ZO 

...EHUM... TABLEAU DE DIMENSION N CONTENANT LES VALEURS DU PROFIL D'HUMIDITE 

* • - ■ I . I 

DO 100 1=1 ,N , ■ ' 

0 EHUM(I)=EO+.66*(TETA(I)-TETO)/BOWEN , ■ ’ . 

MX=3 . ' ■ 

RETURN • ! : 

END 
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SUBROUTINE EOPTI ( A1 , A2, B, NP, EHUM,Z , W,T2W2 , AL, I J) 

DIMENSION EHUMd ),Z(1 ),W(1 ),T2W2(1 ),X(15) • 

EN ENTREE A1 EST UNE CONSTANTE. 

A2 EST UN PARAMETRE CONSTANT LORS D*UN APPEL DU SOUS-PROGRAMME 
B EST UN PARAMETRE CONSTANT EGAL AU RAPPORT DE BOWEN ■ 

NP EST LE NOMBRE DE POINTS DU PROFIL 
EHUM DIMENSION NP CONTIENT PROFIL D'HUMIDITE 

Z DE HAUTEUR 

EN SORTIE W omega (CF.TEXTE) 

T2W2 ...LA DERIVEE SECONDE DE 

TETA PAR RAPPORT A OMEGA) 

R=0.66*A2/B 
DO 300 1=1 ,NP 

IF(IJ.EQ.2)G0 TO 110 : ■ . ' . 

I IF(AL)110,110,120 ■ ' , 

;i CL=-0.02*AL 

' IF(ABS{Z(I)-CL) .LE. 0.000001 )Z(I)=CL 
IF (IJ.EQ.2)G0 TO 111 

IF(Z(I)-CL) 111 ,111 »112 , 

' W(I)=Z(I)*(EHUMn)-R)/A1 ' . 

T2W2 ( I )=-A2* ( 1 .+R/EHUM( I ) )*{ A1 /Z ( I )/EHUM( I ) )**2 
GO TO 300 
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2 W(I)=R*Z(I)**.75*(1 .33*AL/Z(I)-.34)/(2.7*(AB5(AL))**.75) 

W(I)=EHUM(I)-W(I) 

W(I)=W(I)*Z(I)/A1 

' X(I)=R*(1 .+.6»Z(I)/AL)/(EHUM(I)*2.7*(Z{I)/ABS(AL))**.25)+.75 

X(I) = .2*(Z(I)/ABS(AL))*».75+(1 .+. 6^^Z{ I )/AL*X( I ) )/( 2 . 7* (Z { I)/AB5 ( AU)** .25 ) 

' T2W2(I)= X(I)*(-A2*(A1/Z{I)/EHUM(I))**2) . 

. GO TO 300 

,0 DL=0.Q01*AL 

! IF{AB5(Z(I)-DL) .LE. 0.000001 )Z(I)=DL 

IF(Z(I)-DL)111 ,111 ,122 . , 

2 W(I)=R*2*((Z(I)/AL)**.1 )*(.9+.3*Z(I)/AL) 

! W(I)=EHUM(I)_W(I) ' , 

W(I)=Z(I)*W(I)/A1 

, X(I)=R*2*(1 . + .6*Z(I)/AL)*(Z(I)**.1 )/EHUM(I)*(AL**.1 )+1 .1 

, X(I)=1 .2*(Z(I)/AL)**1 .1+(1 .+6*Z(I)/AL)*2*({Z(I)/AL)**.1)*X(I) 

T2W2(I)= X(I)*(-A2*(A1/Z(I)/EHUM(I))**2) 

;0 CONTINUE ' 

RETURN . 

^ END . , ' ' 


89 



APPENDIX II 


TABLE OF NUMERICAL VALUES OP 




(mW/c 



U-J Q — 15 m/s 


: L (m) 

- 5 

- 16 

- 100 

■ - 500 

- 1000 

+ 1000 

+ 500 

+ 100 

+ 16 

+ 5 

1*^10 -^0 



- 9,17 

- 1,46 

- 0,73 

+ 0,75 

+ 1,53 

+ 8,20 



:5o(inW/cm^ 



+ 25,2872 

+ 4,9924 

+ 2,49 

- 2,4841 

- 4,9600 

- 24,4769. 

> ■ • 

' 

: Z{m) 

; 10-3 



+ 1,09 

+2,85.10“^ 

+ 1 ,51 .10“* 

- 1,77.10“* 

- 3,61 .10“*-- 

- 1,88 



• 10-^ 


- 

+ 1,68 

+ 1 ,03.10“ 

+ 4,55.10“^ 

- 6,56.10“^ 

- 2,01 .10“^ 

- 1,91.10“* 

■ 


: 10-^ 



- 6,62.10"^ 

- 1,39.10“'* 

- 7,19.10“^ 

•+ 7,30.10“^ 

+ 1 ,40.10“* 

+ 6,32.10“* 

- 


: 1. 



- 2,66.10“'* 

- 5,91 .10"^ 

- 2,96.10“^ 

+ 2,79.10“^ 

+ 6,01 .10“^ 

+ 3,20.10“* 



: 5 

- 


- 9,02.10“^ 

- 1 , 92 . 1 0“^ 

- 9,35.10"^ 

+ 9,37.10“^ 

+ 1,97.10”^ 

+ 1,09.10"* 

■- 

' - 

: 10 


• ■ - 

- 1,14.10“^ 

- 2,28.10“^ 

- 1,12.10“^ 

+ 1,38.10“^ 

+ 2,96.10“^ 

+ 1 ,79.10“* 




VO 

M 




















IV) 


= 12 m/8 


Urn) 

- 5 

- 16 

- 100 

- 500 

- 1000 

+ 1000 

+ 500 

+ 100 

+ 16 

+ 5 

010 -®o 

% 


- 5,41 

- 0,87 

- 0,44 

+ 0,44 

+ 0,91 

+ 4,04 



5^(mW/cm^ 



+ 10,4026 

+ 2,0556 

+ 1,0263 

- 1,0232 

- 2,0432 

- 10,0^22 

1 .. r , 













ZTm7 

10-3 

c. 


+ 5,54.10“'' 

+ 1,28.10“^ 

+ 6,82.10“2 

- 7,62.10“2 

- 1,56.10“^ 

- 8,27.10“^ 


- 

10-2 



- 1,43.10“'' 

- 3,18.10“2 

- 1,69.10“2 

+ 1 ,64.10“2 

+ 2,9D.10“2 

+ 9,32.10“2 

■ ' •• 


; 10-^ 


• ' 
• 

- 4,78.10“'' 

- 9,83.1D“2 

J 

1- 5,04.1 0“2 

+ 5,02.1D“2 

+ 9,76.10“2 

+ 4,48.10“^ 

' ■ 

. ■ • - 

: 1 


• 

- 1 ,78.10“'' 

- 3,68.10“2 

- 1,84.10“2 

+ 1,74.10“2 

+ 3,72.10“2 

+ .1 ,97.10“'' 



\ 

: 5 



- 6,03.10“2 

- 1 ,12.10“2 

- 5,54.10“^ 

+ 5,73.10“^ 

-2 

+ 1,20.10 

+ 6,58.10“2 

' 


: 10 

- 


- 6,2B.10“2 

- 1 ,32.10“2 

- 6,50.10“^ 

+ 7,97.10“' 

+ 1,70.10"2 

+ 1,02.10“^ 






= B m/8 


L (m) 

- 5 

- 16 

- 100’ 

- 500 

- 1000 

+ 1000 . 

. + 500 

+ 100 

+ 16 ; 

© 1 0 - ©0 

• 

- 11,43 

^ 2,01 

- 0,34 

- 0,17 

+ 0,17 

+ 0,34 

+ 1,82 

: 

+ 12,04 

5jj(mW/cm^ 


+ 13,0612 

+ 1,9578 

! + 0,3075 

+ 0,1935 

- 0,1930 

- 0,3855 

- 1,9074 

-11,100 

Z(mT 

1D-3 


6,17.10"^ 

+ 1 ,27.10"'’: 

+ 2,85.10"^ 

:+ 1,47.10"^ 

- 1 ,55.10"^ 

- 3,23.10"^ 

_1 ,78.10"’ 

1,52 • 

• « 

1 

1 

O 1 

1 1 

ro 1 

1 


7,01 .10"^ 

- 1, 20.10"’ 

- 2,79.10"^ 

- 1,45.10"^ 

1 ,40.10"^ ■ 

4- 2,69.10"^ 

1,17.10"’ 

2,93.10"’ 

; ,0-'. : 


1,05 

- 2,00.10"’ 

- 4,17.10"^ 

- 2,12.10^^ 

2,00.10"^ 

4,07.10"^ 

1,90.10"’ 

+ 1,14 

; 1 


- 4,49.10"** : 

- 6,04.10"^ 

- 1,36.10"^ 

- 6,92.10"^: 

6,60.10"^ 

< 

1 ,39.10"^ : 

7,36.10"^ 

4,54.10"’ 

5 


1,26.10"^ 

2,13.10"^ 

-3 

4,09.10 

' -3 

- 2,02.10 , 

+2,12.10"^ 

4,42.10"^ 

+ 2,40.10"^ 

-+ 1,70.10"’ 

: 1C 


- 1,07.10"^ 

2 *- 4,47*10"^ 
- 1 ,99.10 . 

i * 

- 2,22.10"^ 

-3 

+2,73.10 

w- 5,82.10"^ 

w- 3,49.10"^ 

^ 3,37.10"’ 


+ 5 : 


VO 

VjO 






























- 13,51 


5^(mW/cm^ * 9,8031 


.16^ : 

- 100 : 

- 500 : 

- 1000 : 

+1000 : 

+ 500 : 

+ 100 

4,19 : 

: 

- 0,73 : 

- 0,13 : 

- 0,07 : 

: 

+ 0 1 0 6 • • 

+ 0,12 : 

+ 0,66 


+ 2,6630 


1 . 00.10 


+ 0,4000 


+ 2,36.10 


+ 0,0794 


+ 5,35.10 


+ 0,0396 


- 0,0395 


-3' _3 

+ 2,71.10 •_ 2,90.10 


- 0,0790 


- 6,25.10" 


0,3911 


_ 3,70.10” 


+ 4,40 
- 2,2849 


- 3,90.10 : 


4 50.10“^- 0,35.10 :-1, 71. 10" 8,59.10" 0, 37.10" 1,64.10 :+ 7,75.10 


5,53.10"^- 1,71.10":- 2,67.10 5,26.10"^- 2,64.10" *.. 2,53.10-":- 5,31.10 :+ 2,78.10 


1,51 .10’^]- 7,68.10"'^*.+ 8,02.10"'^:- 1,66.10"^:+ 0,91.10-":+ 6,35.10 ^:+ 2,12.10 


L 1,24.10"'':- 3,87.10"^:- 7,00.10"^:- 1,61.10"^:- 7, 99.1 0”'':+ 9, 81 .1 O"^ :-2 , 1 0 .1 0 :+ 1,25.10 


+ 1,19.10 









2,8 m/s 


“10 “ 


: L (m) 

- 5 

- 16 

• 

• 

: - 100 ■ 

- 500 

• 

- 1000 

+ 1000 

+ 500 

+ 100 

+ 16 

+ 5 

= ®1D 

. - 3,70 

- 1,15 

» 

- 0,20 

- 0,04 

- 0,02 

+ 0,02 

+ 0,03 

+ 0,18 

+ 1,20 

+ 3,67 

S^CmW/cm? 

■ 

1 ,1572 

0,3198 

' 0,0485 

0,0096 

0,0048 

- 0,0048 

- 0,0096 

- 0,0474 

- 0,2793 

- 0,7525 

Z(m) 

10-3 

- 6, 62.10-2 

- 1,75.10“2 

- 1,95.10-3 

4 

- 3, 37. 10"^, 

- 1,72.10“'' 

+r,t3o.jo“^ 

+ 9,12.10"^ 

- 9,7.10“" 

- 4,51 .10“^ 

- 3,12.10"! 

10-2 

' -1 
- 4,59.10 

_1 

- 1,52.10 

* 

- 2,55.10“2;‘_ 5,22.10“^ 

- 2,62.10“^ 

+2,49.10“^ 

+ 4,88.10“^ 

+ 2,3.10“2 

+ 1 ,02.10*"' 

+ 1,73.10“'; 

IQ-^ 

- 4,72.1Q-^ 

- 1 ,45.10“' 

-2^ 

- 2,50.10 

-3 

- 4,99.10 

- 2,48.10"^ 

+2,39.10“^ 

+ 4,73.10“^ 

+ 2,24.10“2 

+ 1,36.10“' 

-1 ' 

+ 3,25.10 . 

1 

- 1,47.10“*' 

- 4,56."0“2 

* -3 

- 7,39.10 

- 1 ,45.10"^ 

- 7,18.10-^ 

+6,89.10“'' 

+ 1 ,45.10“^ 

+ 7,55.10“^ 

+ 4,68.10“^ 

+ 1,0B.10“': 

5 

- 4,17.10“2 

- 1,24,10“2 

- 2,04.10“^ 

- 4,12.10“^ 

- 2,06.10"'* 

+2,16.10“^ 

+ 4, 44.10“" 

+ 2,36.10“^ 

+ 1 ,67.10“^ 

-2 ■ 

+5,60.10 

10 

- 3,11 .10“2 

- 9,8.10“^ 

- 1 ,76.10“^ 

- 4,23.10“^ 

- 2,08.10“^ 

• 

• 

+2,53.10“'' 

• 

• 

+ 5, 34.10“" 

+ 3,17.10“^ 

+ 2,99.10“^ 

+ 1,63.10“' 

r 


VD 

U1 



































Figure 2. Variation of U as a function of X and 
a (according to Hldy and Plate 1967). 
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Figure 4. Variation of z . as a function of U 55 (according 
to Shemdln 1967). 
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Figure 6 . Table summarizing the roughness characteristics of the ocean 
surface (according to Wu 1969). 






























Figure 7. Variation of as a function of U 
the ocean (according to Wu 1969)- 






Figure 8 . Comparison of theoretical expressions of Mo.nln-.QbuKhov 
Swlnbank and KEYPS with experimental data, ln~the "case of Instability 
(according to Bernstein I966). 
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Figure 10. Region of variation of the experimental determinations of the eralsslvlty 
function (work of Brooks, Robinson, Brulnenberg, Elsasser and Kuhn). 





situations : inversion and over-adiabatic conditions 


(u^g ■ 2,B m.«, O o ■ 20* C) 


106 




Figure 12. Characteristic results of the numerical calculations 

obtained for A « 2Q0 C - u.„ • 2,0 m.ml and different stability 

o 1 u _j 

regions . 
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Figure l4. Analytical representation of the first 
derivative 
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Figure 17 . Results obtained for d z j using the simplified 
logarithmic profile expression. Comparison between numerlca: 
calculation and analytical development ( 0 _ ■ 20* C - u^q • 
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Figure 18 . Results obtained for $ using the slm- 

pllfled logarithmic profile expression. Comparison, between 
numerical calculation and analytical development 

(do* 20* C - U1Q - 2,0 m.a) j 
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Figure 19 . Calculated variations of the turbulent flux S between 0 and 10 m 

as a function of the friction velocity of u* (for various humidities and stability 

regions and the emissivities of Brooks) 
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APFENDIX^ I I 


NUMERIC.AL- CALCULATION., PROGRAM (FORTRAN. IV) 


DIMENSION U(5),Z0M(5),DBUKG(1Q) 

DIMENSION Z(15),X(15),TETA(15),EHUM(15),W(15),T2W2(15),\/(15),Y(15) 
DIMENSION WR(10),ZR(10),DIV(10),FI(10),CUMUL{10),T3W3(15),EMIS{15) 

FORMAT (////, 3X-, WIND.. VELO.CITY, AT 10 M.., U10= »,F9..1,' 

^ r. . . ZD= ' ,'E 9 . 3 , * . . L . MffN IW-ORGKHOV^ITTTTTTrtl^'VF 970) — 

FORMAT(/ ,1H , ' .USTAR= • * , F1 0 .4 , • , FLUX QH=» ,F10.4, ' ,EVAPO= 

1 ' ,F10.5, * ,BOWEN= ',F10.4,» ,A= ,F10.5) 

FORMATdH ,,100('*' )) 

^ ' height.- . 9- temperature * 

1 * w ■ • - *- ■ "dZT/DWZ *•) . 

FORMATdH ,5(E15.8,4X,1H*) ) 

F0RMAT(1 HO, 'NP-R * WR * ZR * INTEGRAL' * 

^DIVERGENCE * CUMUL *') 

FORMATdH, I4,1X,»*',5(E13.7,1X,'*’)) 


Initialization bE9;c.0N3TANT 9 ana^ ^^ VALtiks^' 

N=1 4 
AK=0,4 

TET0=293. , 

U(1)=2.8 

U(2)=5. ■ ■ 

U(3)=8. 

U(4)=12. 

U(5)=15. ' ' . 

Z0M(1 )=0. 00001 . 

Z0M(2)=0. 00005 

Z0M(3)=0.0001 

Z0M(4j=0.0005 

Z0M(5)=0.001 . 

DBUKOd )=-1000. 

DBUK0(2)=-500. • , 

]BUK0(3)=-100. . / 

3BUK0(4)=-16 

]BUK0(5)=-5 

)0 100 1 = 6,10 

■=11-I 


)BUKO(I)=-OBUKO(J) 

:0=23.27 

.1=13000 

MX IS THEIJUMBER oi^^THE .-output. UNIT I 



ROUGHNESS, 


_HyMXD.rTY^''I 




STANDARD HEI GHT VAL UES I 
X(1)=Q. • 

X(2)=0. 00001 

X(3)=0. 00005 ' 

X(4)=0.0001 

X(5)=0.0005 

X{6)=0.001 

X(7)=0.005 

X{8)=0.01 

X(9)=0.05 > , , 

X(10)=0.1 

X(11)=0.5 ' 

X(12)=1. 

X(13)=5. 

X(14)=10. ... 

,X(15)=100. 

VALUES OP : THE; .EMiS.SI¥,I.T.Y-.CIJRp:_.PA^^ 

. 80=0.7875 
81=0.1420 
82=0.0065 

UP "to here the NUIVIBER ’ of. P.ATR_^JUl_d ,Z0'') _STUpiED, USED POR^EVALUAT 
’D0 130i=1,5 ^ • 

U10=U(I) . 

Z0=Z0M(I) 

WE CAN SEE THAT THE NUMBER OF PROFILE POINTS AND THE 
POSSIBLE HEIGHT VALUES-,- NEVER ZER0 _ ' 

NP=N+1-I 
DO 120 J=1 ,NP 

K=I+J ■ ■ . ■ ‘ 

20 Z(J)=X(K) 


CALCULATIONS FOR FIXED VALUES OF THE. .MONIN-OBUKHOV 
LENGTHS ARE' BEING .CARRIED. ..O.UT. 


DO 130 J=1 ,10 

IH=1 ^ 

DO 130 IJ=1 ,2 , 

AL=0BUK0(J) 

WRITE(3,1 )U10,Z0,AL 
UET0I=U5TAR(AK,U10,AL,Z0, INDIO) 

• "calculation OF' :THE.- TEMPEIRAT-URE . PROFILE . ] . 


CALL protp(al,ak,uetoi,zo,teto,np,z,teta,ix,ij) 

TET10=TETA(NP) 


CALCULATION OF THE BOWEN RATIO 


B=B0WEN(TET10,TET0) 




-CALCULATION OF THE HEAT’ FLUX :AND. EVAPORATION 
CALCULATION..-OE „THE., HUMIDI-TY-^PROEILE— . - 


CALL FLUEV(UETOI,AL,B.,QHQ,EVAPO) 
CALL PRQHU(EO,B,NP,TETA,TETO,EHUM) 


C ALCULATION. OF .ELE-FEN.TS, INyQLVEB.JD ljriIE.,.. RABIA TXGM - F LUX 


DIVERGENCK_ 


A2=TET0*UET0I*UET0I/(AK*AK*AL*9.81 ) 

CALCULATION OF THE TWO TABLES CONTAINING W AND THE SECOND 
DER IV AT I VE.-QJ’.. TH ET. A„W .IT H- RESPEGT - H.O-.-W- ^ 


CALL E0PTKA1 , A2 , B, NP,EHUM, Z, W, T2W2 , AL, IJ) 

WRITE(MX,2)UET0I,QH0,EVAP0,B,A2 

WRITE(MX,3) 

WRITE(MX,4) 

WRITE(MX,3) 

DO 135 K=1,NP 

35 WRITE(MX,5)Z(K),TETA(K),EHUM(K),W(K),T2W2(K) ■ 

^ -I 

' CALCULATION _OF THE. VALUES ANp.-_WR_ j 

DO 140 K=1 ,10 
L=(MP-K+1 
M=11-K 
WR(M)=W(L) 
iO ZR(M)=Z(L) 

CALCULATION OF EMISSIVITIEs| 


t)8 DO 1 50 M=1 , 1 0 
DO 1 63 K=1 ,MP 
V(K)=W(K)-WR(M) 

IF(V(K))164,162,164 
52 EMI5(K)=0. 

GO TO 163 

^4 IF (ABS(V(K) )-1 .E-4) 165,165,161 

b IF(IH.NE.3) GO TO 300 ■ . . ^ . 

IX=2 

EMI5(K)=0,21 B*(1 .+1 .125*(50RT(AB5(V(K) ) )-7.9E-2)/(50RT(AB5(V(K) ) )+1 7.9E-2) ) 
ED TO 163 

85 IF{IH.i\lE.3) GO TO 301 
IX=2 

EMI5{K)=0.260*AL0G (1 .+1 .043E3*(AB5(V(K) ) ) ) 

GO TO 163 

111 EMIS{K)=0.O55*ALOG(1 .+7.E3*AB5(V{K))) 
h T3W3(K)=T2W2(K)*(EMI5(1 )-EMI5(K)) 


INTEGRATION OF THE FUNCTION T3W3 TABULATED IN THE INTERVAL V(l) TO V(NP)J 
THE INTEGRAL EQUALS FI.SSP PAGE 289 

CALL QTFG(V,T3W3,Y,NP) ^ 

FI(M)=Y(NP) / 



C A LCUL ATIO N OF -T H E D IVERGENCE AT LE VEL jR_] 

L=MP-M+1 

0 DIV{M)=-22.4E-9*TET0**3*EHUM(L)*FI(M)/A1 

' CALCULATION OF -CUMU LATIVE' AREAS. OF DIV OF TH E FUNCT^ON_ OP_^ 
CUMULd )=1 .E-3*DIV(1 ) 

DO 170 M=2,1Q ' . 

CUMUL(M)=CUMUL(M-1 ) + (ZR(M)-ZR(M-1 ))*(DIV(M)+DIV(M-1 ))*.5 , ' 

WRITE(MX,6> 

DO 180 M=1,10 . ■ , 

VIRITE(MX,7)M,WR(M),ZR(M),FI(M) ,DIV(M) ,CUMUL(M) 



; SUB- PRO GRAMS ^ 


FUNCTION U5TAR(AK,U10,AL,Z0, INDIE) 


CALCULATION 'OF 'THE"FRICTION' VELOCITY USTAR(U STAR) 

. . . AK . . . KARMAN CONSTANT 
1 . . .U10...WIND VELOCITY AT 10 METERS 
. . .AL. . .MONIN-OBUKHOV LENGTH 
. . . ZO . . . ROUGHNESS PARAMETER 

I..INDIC. .0 IP THE CALCULATION IS ■ NUMERICALLY POSSIBLE 

..1 IP THE TERM IN THE LOGARITHM IS ZERO, THIS POINT COULD BE 
IMPROVED BY. SPECIFYING THAT THE TERM HAVE AN ORDER HIGHER 
THAN THAT OP EPS 


I'iDIC=0 

iPS=0. 

<‘^L*(EXP(10./AL)-1 .)/Z0 
!^(X-EP5)100,100,1 10 

‘t'<31C»1 

^nURN 

■^fARoAK*U10/AL0G(X) 

:t?URN 

■I'D 





SUBROUTINE PROTP ( AL, AK, USTAR , ZO, TETO, N,Z , TETA , IX, IJ ) 
dimension Z(1 ) ,TETA(1 ) 


OATHITATTON OF A TEMPERATURE PROFILE THETA(Z) IS A FUNCTION OF Z 
5mLE CM?AimNG ™ VARIOUS VALUES OF THE OBSERVATION HEIGHTS, 

DIMENSION N 

..AL MONIN-OBUKHOV LENGTH 

..AK KARMAN CONSTANT 

. .USTAR. FRICTION VELOCITY 
ZO ROUGHNESS PARAMETER 

■ M numW of different height values 

'.THETA '.TABLE CONTAINING THE RESULTS OF DIMENSION N 
. .TETO. .FIXED PROFILE VAL^UE_POR ^0 - 


GO TO ( 201, 202), I J 
11 T=ABS{AL*EXP(-0.02)-1 .) 

■ IX=1 
G=9.B1 

T=TET0*(1 .+U5TAR*U5TAR*AL0G(T/Z0)/AK/AK/G/AL) 

DO 90 1=1 ,N 
IF(AL)100,100,121 
iO BL=-.02*AL 

IF(AB5(Z( I )-BL),EL. 0.000001 )Z{I)=BL ■ 

IF(Z{I)-BL)110,11C,120 

D CL=AB5(BL) . , , 

TETAU )=T+,2.8*U5TAR*USTAR*(1 ./Z(I)**.5-1 ./CL** .5)/G/AB5(AL)** .5 

GO TO 90 

0 5=AB5{AL*(EXP(Z(I)/AL)-1.)) 

IF (ABS(5) .L£. 0.00001 )TETA(I)=TETO 
IF(TETA(I) .EQ.TETO)GO TO 90 

TETA( I )=TET0*(1 .+U5TAR^^U5TAR*AL0G(5/Z0)/AK/AK/G/AL) 

GO TO 90 

1 DL=0.001*AL 

IF(ABS(Z(I)-DL).LE. 0.000001 )Z(I)=DL 
IF(Z(I)-DL) 110,110,122 

2 TPRIM=AB5{AL*(EXP(0.001 )-1 ) ) 

jPrqM=TETO* ( 1 . +U5TAR*U5TAR*AL0G ( TPR IM/ZO ) /A K/AK/G/AL ) 

TETA( I )=TPR0M220.*USTAR*U5TAR*TET0*( (Z ( I )/AL)**0 .1 -0.001 1*’*0 . 1 )/AK/AK/G/AL 
;0 CONTINUE 
,0 RETURN 

2 TETA(I)=TET0*(1 .+U5TAR**2*AL0G(Z(I)/Z0)/AK/AL/AK/G) 

IX=2 , , 

RETURN 
: END 



FUNCTION B0WEN(TET_1 0,TET0) _ _ 

CALCULATION OF THE BOWEN RATIO 

DELTA=TET10-TET0 
IF(DELTA)100,100,110 
0 B0WEN=.1-.09*DELTA 

■ mx =3 'n 

. RETURN / ^ 

n BnWFN- .R*nFLTA 


SUBROUTINE FLUEV( U5 TAR , AL, BOWEN, QHO, EVAPO) 


THIS SUBROUTINE CALCULATES 

. . QHO FLUX OP SENSIBLE HEAT 

. .EVAPO. EVAPORATION 
WHEN THE FOLLOWING ARE- KNOWN 
. .USTAR. FRICTION VELOCITY 

. .AL NON ZERO MONIN-OBUKHOV LENGTH 

..BOWEN. NON ZERO BOWEN. RA.TIO . 


QHO=- . 9E+Q4*U5TAR*U5TAR*USTAR/AU 

EVAPp=1 .5E-02*QH0/B0WEN 

MX=3' 

RETURN 

END 


SUBROUTINE PROHU (EO, BOWEN, N, TETA , TETO, EHUM) 

DIMENSION TETA(1 ),EHUM(1 ) 

' CALCULATION 0? THE.. HUMIDIT.Y.. PROFILE _JF THE .TEMPERATURE PROFILE IS KNOWN 
■ . . .EO. . . . .yALUE_-O.E^THE_.HUMini.TY--F-OE-Z-=Z.O. ROUGHNESS PARAMETER 

. BOWEN:.....BOWENvJIA.TXO_. . 

...-■.-N,.. . .■.-.NUMBER OF- DIFFERENT' HEIGHT VALUES 

. . .. r'TSTJlj-. N DIMENSION TABLE CONTAINING THE TEMPERATURE PROFILE VALUES 

.... .TETO. . .CONSTANT CORRESPONDING TO THE VALUE OP ZO 

. . .EHUM.'.VN DIMENSIONAL TABLE CONTAINING THE HUMIDITY PROFILE VALUES 


DO 100 1=1 ,N 

EHUM( I ) =E0+ . 66* ( TETA { I ) -TETO ) /BOWEN 

HX=3 

RETURN 

END 


SUBROUTINE E0PTI(A1 , A2 , B, NP, EHUM, Z , W, T2W2, AL, I J) 

DIMENSION EHUM(1 ),Z(1 ),W(1 ),T2W2(1 ),X(15) ' 

\ . 


AS INPUT 


AS OUTPUT 


A1 IS A CONSTANT. 

A2 IS A CONSTANT PARAMETER WHEN THE 'SUB-PROGRAM IS CALLED 
B IS A CONSTANT PARAMETER EQUAL TO THE BOWEN RATIO 
NP IS THE NUMBER OF PROFILE POINTS 
EHUM DIMENSION NP CONTAINS THE HUMIDITY PROFILE 
Z HEIGHT 

W OMEGA (SEE TEXT) 

T2W2 SECOND DERIVATIVE OF 

. ... TE.TA -WITH-RESPECX TO OMEQA) 


»'0.66*A2/B 
jO 300 1=1 ,NP 
!niJ,EQ.2)GO TO 110 
I(AL)110,110,120 
‘I=-0.02*AL 




SHAD5(2(I)-CL) .LE. 0.000001 )Z(I)=CL 
(IJ.EQ.^0 TO 111 
;tfU(I)-CL) 111,111,112 



2 W(I):=R*Z(I)»*.75*(1 .33*AL/Z(I)-.34)/(2.7*(ABS(AL) )**.75) 

W(I)«EhUM(I)-W(I) 

W(I)=W(I)*Z(I)/A1 

X(I)=R*(1 .+.6*Z(I)/AL)/(EHUM(I)*2.7*(Z(I)/ABS{AL))**.25)+.75 

X(I) = .2*{Z(I)/AB5(AL) )*^^.75+(1 . + . 6*Z ( I )/AL*X( I ) )/{2 . 7* (Z ( I )/AB5 ( Aa)*» .25) 

T2W2(I)= X(I)*(-A2*(A1/Z(I)/EHUM(I))**2) 

GO TO 300 
|DL=D.001>AL 

IF(AB5(Z(I)-DL) .LE. 0.000001 )Z(I)=DL 
IF(Z(I)-DL)111,111,122 
^W(I)=R*2*((Z(I)/AL)**.1)*(.9+.3*Z(I)/AL) 

W(I)=EHUM(I)-W(I) 

W(I)=Z(I)^W(I)/A1 

X{I)=R*2*(1 . + .6*Z(I)/AL)*(Z(I)**.1 )/EHUM(I)*(AL**.1 )+1 .1 
X(I)=1 .2*(Z(I)/AL)**1 .1+(1 .+6*Z(I)/AL)*2*((Z(I)/AL)**.1 )*X(I) 

T2W2(I)= X(I)*(-A2*(A1/Z(I)/EHUM{I))**2) 

CONTINUE 

RETURN ^ 

END 
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